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Abstract 

We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to 
produce six temperature maps from the first 473 days of Planck-RFl survey data. These maps provide an accurate rendition of the sky emission at 
I 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9!! to 4!6. The detector noise per (eff'ective) beam solid angle is 
. • respectively, 10, 6 , 12 and 39 yt/K in HFI four lowest frequency channel (100-353 GHz) and 13 and 14 kJy sr~^ for the 545 and 857 GHz channels. 
^ Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5 % and the 353 GHz relative calibration 
is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB 
rS power spectrum (50 < ^ < 2500), are intercalibrated at better than 0.2 %. 

Key words. Cosmology: cosmic background radiation - Surveys - Methods: data analysis 
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1. Introduction 

This paper, one of a set associated with the 2013 release of 
data from the Planck^ mission (Planck Collaboration I 2013), 
describes the processing of data from the High Frequency 



* Corresponding author: F. R. Bouchet, bouchet@iap . f r. 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 



Instrument (HFI) from satellite data to calibrated and charac- 
terised maps. HFI (Lamarre et al. 2010; Planck HFI Core Team 
2011a) observes in the 100, 143, 217, 353, 545, and 857 GHz 
bands with bolometers cooled to O.IK. The HFI instrument 
comprises 50 signal bolometers, as well as two dark bolomet- 
ers, 16 thermometers, a resistor, and a capacitor used for mon- 
itoring and housekeeping. The count of 50 bolometers includes 



telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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twelve polarization sensitive bolometer (PSB) pairs, four each 
at 100-353 GHz; the rest are unpolarized spider-web bolometers 
(SWBs). We describe the steps taken by the HFI Data Processing 
Centre (hereafter DPC) to transform the packets sent by the 
satellite into sky maps at HFI frequencies, with the help of ancil- 
lary data, for example, from ground calibration. These are tem- 
perature maps alone, as obtained from the beginning of the first 
light survey on 13 August 2009, to the end of the nominal mis- 
sion on 27 November 2010. 

Planck defines a sky survey as the time over which the spin 
axis rotates by 180°, a period close to six months in duration 
over which about 95 % of the sky is covered at each frequency. 
During routine operations, Planck scans the sky by spinning in 
circles with an angular radius roughly 85°. The spin axis fol- 
lows a cycloidal path on the sky by periodic step-wise displace- 
ments of 2\ resulting in typically forty (35 to 70) circles of typ- 
ical duration 46 minutes, constituting a stable pointing period 
between repointings. The scanning strategy is discussed in more 
detail in Tauber et al. (2010), Planck Collaboration I (201 1), and 
Planck Collaboration I (2013). The 15.5 months of nominal mis- 
sion survey data then provide 2.5 sky surveys, and maps are 
provided for the first two sky surveys separately as well as for 
the complete nominal mission. As a means to estimate aspects 
of the noise distribution, we also deliver "half-ring" maps made 
out of the first and second half of each stable pointing period. 
Maps are produced for individual detectors, as well as for aver- 
ages over each band and for selected detector sets defined within 
each band (see Table 3). 

The next section provides an overview of HFI data pro- 
cessing. Section 3 is devoted to the processing of Time Ordered 
Information (hereafter TOI) from individual detectors to produce 
cleaned timelines. These timelines are used to estimate the tem- 
poral noise properties in Sect. 3.10 and to determine the detector 
pointings and beams in Sects. 4 and 5. Section 6 discusses the 
creation of maps and their photometric calibration, while Sect. 7 
presents tests applied to assess the consistency and accuracy of 
the products. For completeness, component separation and fur- 
ther processing is briefly described in Sect. 7.5. Section 8 con- 
cludes with a summary of the characteristics of the HFI data 
delivered, as currently processed. 

Some of the specific processing steps of HFI data are 
described more fully elsewhere: Planck Collaboration VII 
(2013) discusses the transfer function and beams; 
Planck CoUaboration VIII (2013) the calibration of HFI 
detectors; Planck Collaboration IX (2013) the determination of 
the spectral bands for each detector and their combination; and 
Planck Collaboration X (2013) the eff'ect of so-called "glitches" 
such as cosmic-ray hits on the detectors. The processing of 
data from the Low Frequency Instrument (LFI) is discussed in 
Planck Collaboration II (2013). Technical details applying to 
specific data products are discussed in Planck Collaboration ES 
(2013). We have applied to the data products delivered many 
consistency and validation tests to assess their quality. While 
the products meet a very high standard, as described here, 
we did find limitations. Their mitigation, and related data 
products, are left to future releases. In particular, HFI analysis 
revealed that nonlinear eff'ects in the on-board analog-to-digital 
converters (ADC) modified the recovered bolometer signal. In 
situ observations over 2012-2013 are measuring this eff'ect, 
and algorithms have been developed to explicitly account for 
it in the data analysis. However, the first-order eff'ect of the 
ADC non-linearity mimics a gain variation in the bolometers, 
which the current release measures and removes as part of the 



calibration procedures. This is discussed further in Sects. 6.2 
and 7.2.1. 

The mapmaking procedure uses the full intensity and polar- 
ization information from the HFI bolometers. The current ana- 
lysis cannot guarantee that the large-scale polarization signal is 
free from systematic eff'ects. However, the preliminary analysis 
shows that the small-scale maps have the expected CMB con- 
tent at high signal to noise, as discussed in Sect. 6.7 below, and 
in Planck Collaboration I (2013) and Planck Collaboration XVI 
(2013). Although we do not use these maps for cosmological 
measurements, future work will use them for investigations of 
the properties of polarized emission of the Galaxy. 



2. HFI data processing overview 

The processing of HFI data proceeds according to a series of 
levels, shown schematically in Figure 1. Level 1 (LI) creates a 
database of the raw satellite data as a function of time (TOI ob- 
jects). The full set of TOI comprises the signals from each HFI 
bolometer, ancillary information (e.g., pointing data), and asso- 
ciated housekeeping data (e.g., temperature monitors). Level 2 
(L2), the subject of this paper, uses these data to build a model 
of the HFI instrument, the Instrument Model (IMG), produces 
cleaned, calibrated timelines for each detector, and combines 
these into aggregate products such as maps at each frequency. 
Level 3 (L3) takes these instrument- specific products and de- 
rives various products: component- separation algorithms trans- 
form the maps at each frequency into maps of separate astro- 
physical components; source detection algorithms create cata- 
logues of Galactic and extragalactic objects; finally, a likelihood 
code assess the match between a cosmological and astrophysical 
model and the frequency maps. 

Of course, these processing steps are not done completely 
sequentially: HFI data is processed iteratively. In many ways, 
the IMG is the main internal data product from Planck, and the 
main task of the HFI DPC is its iterative updating. Early ver- 
sions of the IMG were derived from pre-launch data, then from 
the first light survey of the last two weeks of August 2009. 
Further revisions of the IMG, and of the pipelines themselves, 
were then derived after the completion of successive passes 
through the data. These new versions included expanded inform- 
ation about the HFI instrument: for example, the initial IMG con- 
tained only coarse information about shape of the detector angu- 
lar response (i.e., the full-wifth at half-maximum of an approx- 
imate Gaussian); subsequent revisions included full measured 
harmonic- space window function. 

In somewhat more detail, LI software fills the database 
and updates, daily, the various TOI objects. Satellite atti- 
tude data, sampled at 8 Hz during science data acquisition 
and at 4 Hz otherwise, are resampled by interpolation to the 
180.37370 Hz (hereafter 180.4 Hz) acquisition frequency of the 
detectors, corresponding to the integration time for a single 
data sample; further information on LI step were given in 
Planck HFI Core Team (2011b). Raw timelines and housekeep- 
ing data are then processed by L2 to compensate for instrumental 
response and to remove estimates of known artefacts. The vari- 
ous steps in TOI processing are discussed in Sect. 3. First, the 
raw timeline voltages are demodulated, deglitched, corrected for 
the bolometer non-linearity, and for temperature fluctuations of 
the environment using correlations with the signal TOI from 
the two dark bolometers designed as bolometer plate temperat- 
ure monitors. Narrow spectral lines caused by the "^He-JT (4K) 
cooler are also removed before deconvolving the temporal re- 
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Figure 1. Overview of the data flow and main functional tasks of the HFI Data Processing Centre. Level 1 creates a database of 
the raw satellite data as a function of time. Level 2 builds a data model and produces maps sky maps at the six frequency of the 
instrument. This flow diagram illustrates the crucial role of the Instrument Model (I MO), which is both an input and an output of 
many tasks, and is updated iteratively during successive passes of the data. Level 3 takes these instrument-specific products and 
derives final astrophysical products. This paper is mostly concerned with level 2 processing and its validation. 



sponse of the instrument. Finally, various flags are set to mark 
unusable samples. 

Further use of the data requires knowledge of the pointing 
for individual detectors as discussed in Sect. 4. During a single 
stable pointing period, Planck spins around an axis pointing 



towards a fixed direction on the sky (up to an accounted-for 
wobbling), repeatedly scanning approximately the same circle 
(Planck Collaboration I 2013). The Satellite is re-pointed so that 
the spin axis follows the Sun, and the observed circle sweeps 
through the sky at a rate of approximately one degree per day. 
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Assuming a focal plane geometry, i.e., a set of relations between 
the satellite pointing and that of each of the detectors, we build 
rings of data derived by analysing the data acquired by a detector 
during each stable pointing period ("ring" refers to the data 
obtained during a single stable pointing period). This redund- 
ancy permits averaging of the data on rings to reduce instrument 
noise. The resulting estimate of the sky signal can then be sub- 
tracted from the timeline to estimate the temporal noise power 
spectral density, a useful characterization of the detector data 
after TOI processing. This noise may be described as a white 
noise component, dominating at intermediate temporal frequen- 
cies, plus additional low and high-frequency noise. The effect 
on maps of the low frequency part of the noise can be partially 
mitigated by determining an offset for each ring. These so-called 
"destriping" offsets are obtained by requiring that the difference 
between intersecting rings be minimised. Once the offsets are 
removed from each ring, the rings are co-added to produce sky 
maps. 

As explained in Sect. 6, a complication arises from the fact 
that the detector data include both the contribution from the 
Solar dipole induced by the motion of the Solar System through 
the CMB (sometimes referred to as the "cosmological" dipole), 
and the orbital dipole induced by the motion of the satellite 
within the Solar System, which is not constant on the sky and 
must therefore be removed from the rings before creating the 
sky map. The Solar dipole is used as a calibration source at lower 
HFI frequencies, and bright planet fluxes at higher frequencies. 
Since we need this calibration to remove the orbital dipole con- 
tribution to create the maps themselves, the maps and their calib- 
rations are obtained iteratively. The dipoles are computed in the 
non-relativistic approximation. The resulting calibration coeffi- 
cients are also stored in the IMO, which can then be used, for in- 
stance, to express noise spectra in Noise Equivalent Temperature 
(NET). The destriping offsets, once obtained through a global 
solution, are also used to create local maps around planets. As 
described in Sect. 4 these are used to improve our knowledge 
of the focal plane geometry stored in the previous version of 
the IMO and to improve measurements of the "scanning" beam 
(defined as the beam measured from the response to a point 
source of the full optical and electronic system, after the filtering 
done during the TOI processing step, described in Sect. 5). 

The ring and mapmaking stages allow us to generate many 
different maps, e.g., using different sets of detectors, the first or 
second halves of the data in each ring, or from different sky sur- 
veys. Null tests using difference maps of the same sky area ob- 
served at different times, in particular, have proved extremely 
useful in characterising the map residuals, described in Sect. 7. 



3. TOI processing 

In the LI stage of processing [previously described in Section 3 
of Planck HFI Core Team (2011b)], the raw telemetry TOI are 
unfolded into one time series for each bolometer. The signal is 
regularly sampled at 180.4 Hz. We denote as TOI processing 
the transformation of the TOI coming from LI into clean TOI 
objects, which can be used for mapmaking after focal plane 
geometry reconstruction and photometric calibration have been 
done. The general philosophy of the TOI processing is to modify 
the timelines as little as possible, and therefore to flag regions 
contaminated by systematic effects (e.g., cosmic-ray glitches). 
We deal with each bolometer signal separately. Aside from al- 
lowing the possible flagging of known bright sources, the only 
pointing information that is used in the main TOI processing is 



the phase (see Sect. 3.3) so the TOI processing assumes per- 
fect redundancy of the data within a given pointing period. The 
output of TOI processing is not only a set of clean TOI but 
also accompanying qualifying flags and trend parameters used 
internally for detailed statistics. Moreover, all data samples are 
processed while only clean samples will be projected on maps. 
For beam measurement (see Sect. 5), specific processing is per- 
formed on pointing periods that are close to Mars, Jupiter or 
Saturn (see Sect. 3.11). 

A flow diagram is shown in Fig. 2 illustrating the TOI pro- 
cessing steps detailed in the following subsections. Section 4 
of Planck HFI Core Team (2011b) presents the early version of 
the TOI processing. The pipeline has changed sufficiently since 
to warrant a self-contained global description of the TOI pro- 
cessing. We refer the reader to Planck Collaboration ES (2013) 
and the various companion papers mentioned above for more de- 
tail. The changes mostly reflect the improvement in performance 
and in understanding of the underlying effects. The TOI are not 
delivered in the present data release but their processing is an es- 
sential, though hidden, ingredient in the delivered maps. Some 
of the systematic effects arising in the map analysis can only be 
understood by referring to the bolometer timeline behaviour and 
processing. 

It was recently realised that some apparent gain variations, 
spotted comparing identical pointing circles one-year apart, ac- 
tually originate in non-linearities in the bolometer readout sys- 
tem ADCs. Note that the ADC non-linearity is not corrected for 
in the TOI processing, but rather as an equivalent gain variation 
at the mapmaking level (see Sect. 6.2 and Sect. 7.2.1). 

3.1. Input flags 

Strong signal gradients can adversely affect some stages of TOI 
processing. We flag the data expected to have strong gradients 
using only the pointing information, only in intermediate stages 
and not for the mapmaking stage. We also flag data where the 
pointing is known to be unstable. Input flags come from the fol- 
lowing: 

1. A Point-source flag. It is based on the locations of sources in 
the ERCSC. A mask map is generated around each of them 
from which flag TOI objects are created using the pointing 
information. 

2. A Galactic flag. It is based on IRAS maps with a threshold 
that depends on the frequency. 

3. A BigPlanet flag. Any sample that falls within a given dis- 
tance of Mars, Jupiter or Saturn is retained for this flag. 

4. An Unstable pointing flag (see Sect. 4). This accounts for 
depointings between stable pointing periods and other losses 
of pointing integrity. 

Flagged data are either processed separately or ignored, depend- 
ing on the specific analysis. The flags are described in more de- 
tail in Planck Collaboration ES (2013). 

3.2. Demodulation 

The bolometers are AC square-wave modulated to put the ac- 
quisition electronics l/f noise at high temporal frequencies 
(Lamarre et al. 2010). A demodulation step is done as follows. 
First a one-hour running average of the modulated timeline is 
computed, known as the AC offset baseline. This is done by ex- 
cluding data that are masked due to glitches (see below) or by the 
Galactic flag (see above). Once the AC offset baseline is subtrac- 
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Figure 2. Main TOI processing steps. Pipeline modules are represented as filled squares. TOI: Time Ordered Information. ROI: ring 
ordered information (i.e., one piece of information per stable pointing period). Ring: data obtained during a stable pointing period 
(typically 40 satellite revolutions of one minute each). 



ted from the raw timeline, a simple (+, -) demodulation is ap- 
plied. The overall sign of the signal is set to get a positive signal 
on point sources and Galaxy crossings. This AC off'set baseline 
removal is needed in order to correct for the slow drift of the 
zero level of the electronics. Any possible drift of the baseline 
at a timescale smaller than one hour is dealt with at the filter- 
ing stage (see below). The baseline varies very smoothly over 
the mission and fluctuates by less than ten minimum resolution 
units from the middle of 65536 values allowed by the on-board 
ADC, discussed in greater detail in Sect. 7.2.1. 

3.3. Deglitching and gap-filling 

The timelines are afl'ected by obvious ''glitches" (cosmic ray 
hits and other large excursions) at a rate of about one per 
second. This generates a huge Poisson noise if not dealt 
with. Such glitches are detected as a large positive signal 
followed by a roughly exponential tail. There are two ba- 
sic classes of glitches afl'ecting bolometers. The statistical 
and physical understanding of the two populations is given 
by Planck HFI Core Team (2011b), revised and expanded by 
Planck Collaboration X (2013). Here we present the general al- 
gorithm. 

The type of TOI used in the deglitching is slightly difl'erent 
from the one used in the main processing described in the fol- 
lowing sections. Before deglitching, the timeline is demodulated 
and digitally filtered with a three-point (0.25, 0.5, 0.25) mov- 
ing kernel. Linear interpolation is performed on those parts of 
the timeline when the bolometer pointing nears Mars, Jupiter 
or Saturn, as determined by the BigPlanet flag (see Sect. 3.1). 
This step is done to treat the pointing periods containing a planet 
crossing whose large gradients would otherwise be conflated 
with glitches. 

The algorithm also requires an estimate of the sky signal in 
order to assess the magnitude of excursions about the mean. It 
relies on the fact that the signal component of the timeline is 
periodic within a stable pointing period (up to the slight wob- 
bling of the satellite spin axis). We construct a phase-binned ring 
(PBR), a useful estimate of the sky signal obtained by averaging 
the unflagged TOI samples in bins of constant satellite rotation 
phase. The phase of a sample is given by the pointing reconstruc- 



tion pipeline and varies continuously from to In for each scan 
circle (about one minute). The bin size is about the width of a 
single sample, i.e., V.l and the definition of the zero of the phase 
is irrelevant. 

The algorithm treats each pointing period separately and one 
bolometer at a time. The localization of glitches is done with 
a sigma-clipping method applied to the sky-subtracted TOI. A 
template fitting method is used to identify the type of each glitch. 
After masking and subtracting the series of fitted glitches from 
the original timeline, the PBR is then recomputed as the average 
of unflagged samples. Several iterations (generally six) are per- 
formed until the variation of becomes negligible. The spike 
part of each glitch is flagged and the exponential tail, below the 
3.3 (T level, is subtracted using the last iteration of the glitch tem- 
plate fitting process. 

Most of the samples flagged by this process are due to cos- 
mic ray hits. Figure 3 shows the mean evolution per channel 
of the fraction of flagged samples over the full mission. More 
glitch statistics are given by Planck Collaboration X (2013) and 
Planck Collaboration ES (2013). 

All gaps due to flagged samples within a ring of data are 
replaced with an estimate of the signal given by the PBR. For 
samples around planets, we fill the timeline with the values read 
from a frequency map (made by excluding planets) from a previ- 
ous iteration of the data processing at the corresponding pointing 
coordinates. 

3.4. Bolometer non-linearity correction 

If the environment of a bolometer changes, its response will 
change accordingly. The extreme thermal stability of Planck is 
shown by the measured total power level, which is seen as almost 
constant (see the next paragraph). Nevertheless, the timelines 
must be corrected to account for the slightly varying power ab- 
sorbed by the bolometer coming from the sky load and lOOmK 
bolometer plate temperature fluctuations. For that purpose, the 
voltage-to-power conversion step is assumed to be a simple 
second order polynomial applied to the timeline: 

p = f(i.^), (1) 

Sq\ Vol 
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Figure 3. Evolution of the fraction of the flagged data per bolo- 
meter averaged over a channel for the six HFI channels. A run- 
ning average over 3 1 pointing periods (approximately a day) is 
shown. The peak around day 455 comes from the sorption cooler 
exchange, while the peaks after day 500 are due to solar flares. 



where v is the out-of-balance demodulated voltage (i.e., the 
voltage difl'erence between the square-wave compensating input 
voltage and the total bolometer voltage), 5*0 is the responsivity 
(typically of order 10^ V/W) for the fiducial v = voltage. It is 
measured on ground-based I-V curves, but its exact value is un- 
important as it is calibrated out. The parameter vq is determined 
for each bolometer from measurements made during the CPV 
phase. In general the vq values (5 to 200 mV) are equal to or 
larger than expected using a static bolometer model. The main 
reason is that a TOI sample is an average of 40 on-board meas- 
urements (half a modulation period), which includes an elec- 
tronic transient (see e.g., Catalano et al. 2010). The non-linearity 
correction, as measured by v/vq in Eq. 1, is static and accurate 
for a slowly varying signal, in particular for dipole calibration. 
The term v/vq is of the order of 10"^ at most, over half a dipole 
period, the eff'ect being largest in the 100 GHz band. 

During the mission, the drift of the zero level (as meas- 
ured by v/Sq) is at most roughly 7 fW for bolometers up to the 
353 GHz channels, and 50 fW for the 545 and 857 GHz chan- 
nels, such that the non-linearity correction implies a drift in the 
gain of typically 10"^ at most, for all detectors. The correction 
has a relative precision better than 10%, so that the gain un- 
certainty is at the 10"^ level. The eff'ective gain and uncertainty 
is, however, further aff'ected by the off'set in the ADC electronics, 
as discussed in Sect. 7.2.1 below. The dynamical bolometer non- 
linearity on very strong sources (Jupiter and the Galactic centre) 
are not corrected in the nominal TOI. 

In addition, there are limits associated with the saturation of 
the analog-to-digital converters. As discussed in the explanatory 
supplement, Jupiter and the Galactic centre are the only signals 
on the sky strong enough to trigger this eff'ect. 

3.5. Thermal drift decorrelation 

The removal of the common mode due to temperature fluctu- 
ations of the 100 mK cooler stage is based on measured coupling 
coefl&cients between the bolometers and the bolometer plate tem- 
perature. The coefficients were measured during the CPV phase 
(typically 50 pW/K for bolometers up to the 353 GHz chan- 



nel, and 300 pW/K for the 545 and 857 GHz channels). The 
two dark bolometers are used as a proxy for the bolometer plate 
temperature fluctuations, as HFI 100 mK thermometers have 
too many cosmic ray hits to be used. The two dark bolometer 
timelines are deglitched (in the same way as the other bolomet- 
ers) then smoothed with a one minute flat kernel. From each bo- 
lometer timeline, a linear combination of the two smoothed dark 
timelines is subtracted. During this stage, the timelines of all sig- 
nal bolometers are flagged at the samples where any of the dark 
bolometers is flagged for at least 30 seconds. This is done in 
order to suppress the impact of large glitches happening on the 
bolometer plate. We found empirically that this method excludes 
automatically the worse solar flares and in particular the rising 
common mode induced by massive cosmic ray events. Note that 
the total range of the temperature fluctuation of the bolometer 
plate is less than 80yuK, except for some strong solar flares. As 
shown in Fig. 6 of Planck HFI Core Team (201 la), the temper- 
ature fluctuations of the 1.6 K and 4K stages have a negligible 
influence on the TOI noise properties. 

3.6. 4 K cooler line removal 

The "^He-JT cooler (hereafter, 4 K cooler) produces strong, nar- 
row lines in the power spectral density (PSD) of HFI data. 
These arise from electromagnetic interference of the 4 K cooler 
moving parts with the bolometer readout wires, situated in the 
Planck service module. The 4 K cooler main frequency (40 Hz) 
and its harmonics appear as very narrow lines PSD of the sig- 
nal with an intensity of 10 to 100 times above the white noise. 
The oscillations of the moving parts of the 4 K cooler are syn- 
chronously driven by the same computer that records the data 
samples. Hence for a given stable 4 K cooler line, only one fre- 
quency is affected in the PSD, i.e., a narrow line is produced, 
only broadened by the pointing period integration time. Only 
nine individual frequencies are detectable, each of which can be 
traced to some 4 K cooler harmonic unfolded by the AC modu- 
lation. In the signal domain, they are at 10, 20, 30, 40, 50, 60, 
70, 80, and 17 Hz, which are these fractions of the modulation 
frequency: 1/9, 2/9, 3/9, 4/9, 5/9, 6/9, 7/9, 8/9, 5/27. Some 
thermal instability in the service module makes the line proper- 
ties to change on timescales larger than ten minutes. 

We now show the method used to correct the 4 K cooler line 
contamination in the timelines. In order to prepare the data on 
which to compute the Fourier coeflftcients, we build an interme- 
diate set of signal-removed T0I4K objects. We "unroll" the PBR 
by using as a signal estimate the PBR value from the sky loca- 
tion observed at each time. We then subtract this unrolled signal 
from the original TOI objects. Point sources are linearly inter- 
polated in these intermediate T0I4K. Then, we measure the co- 
sine and sine Fourier coeflftcients on the T0I4K for each of the 
nine lines, once per pointing period. We thus neglect variations 
of these components within the duration of a pointing period. 

We can then subtract from the original TOI objects a timeline 
made of the nine reconstructed components. This notch fil- 
ter scheme produces ripples around strong point sources(see 
Planck HFI Core Team 2011b). The origin of this problem lies 
in the relative position of the 4 K cooler line frequencies with 
respect to the harmonics of the spin frequency. Most of the time 
(typically 97 % of the pointing periods), a given 4 K cooler line 
does not overlap with any spin frequency harmonics and thus 
does not affect the signal at all. The spin frequency is very stable 
within a pointing period and varies from one pointing period to 
another by a factor of at most 10""^ . Hence, in so-called resonant 
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pointing periods, a 4 K cooler line overlaps one harmonic of the 
spin frequency. The Fourier coefficients are perturbed by the sig- 
nal and no longer represent the systematic effect. For those point- 
ing periods only, we interpolate the Fourier coefficients from ad- 
jacent non-resonant pointing periods and subtract the systematic 
timeline from the original TOI objects in the same way as the 
other pointing periods. 

While tests with simulations (described in the explanatory 
supplement) demonstrate that any residual 4 K line contamina- 
tion is reduced a less than 3 % of other sources of noise, these 
lines also affect our ability to characterise and remove the ADC 
non-linearity (discussed in Sect. 7.2.1), and are thus still a sub- 
ject of active analysis. 

The efficiency of our removal may be judged by compar- 
ing map power spectra with this step switched on or off. It 
is found that the affected multipoles are at ^ ^ 60 (//I Hz), 
i.e., 600, 1020, 1200, 1800,... When the module is switched 
on, the line residuals amount to no more than 2.5 % of the 
noise level at those multipoles and much less elsewhere (see 
Planck Collaboration ES 2013) 

3. 7. Fourier transform processing 

The bolometer response suffer from time-constant effects, which 
must be corrected. These effects, modelled as a complex trans- 
fer function whose parameters are determined on planet data 
(see Sect. 3.1 1), are deconvolved at the end of TOI processing. 
An additional filter is applied in order to avoid a large increase 
of noise (without much signal) produced by the deconvolution 
in the last 20 Hz near the modulation frequency. We refer the 
reader to a dedicated accompanying paper for a complete de- 
scription (Planck Collaboration VII 2013). Filtering and decon- 
volution are done during the same Fourier and inverse Fourier 
transform stage. Chunks of 2^^ samples are used at a time with 
a standard fast-Fourier transform (FFT) and inverse transform. 
The deconvolved and filtered data are then computed and saved 
for the middle 2^^ samples. The process is continued by shift- 
ing the input samples by 2^^ until the end of the timeline. This 
ensures continuity in the recovered final timelines. 

For the 545 and 857 GHz channels, the filtering is digital and 
consists of a simple local smoothing kernel: 0.25, 0.5, 0.25. For 
these channels, the filtering is applied first; then the deconvolu- 
tion is done with the FFT module. 

3.8. Jump corrections 

We now build another set of intermediate signal-removed TOI 
from the deconvolved data. These consist mostly of noise, ex- 
cept in regions of strong gradients, for example near the Galactic 
plane. Some pointing periods are seen to be affected by a sudden 
jump (either positive or negative) in the signal- removed TOI. 
We therefore correct for that jump by subtracting a piecewise- 
constant template from the timeline, while preserving the mean 
level. The exact jump location not being precisely fitted, the TOI 
are flagged around the recovered position with 100 samples on 
both sides. 

We find on average 17 jumps per day, over all bolomet- 
ers, with jumps affecting a single bolometer at a time. Hence, 
a fraction less than 10"^ of data is lost in the flagging process. 
The jump rate fluctuates during the mission with a peak-to-peak 
value of nine jumps per day. So far, there is no real explanation 
for these events, although we can suspect cosmic ray violent hits 
on the warm electronics. 



We have checked the jump correction process on simulated 
jumps of various intensities added to pointing periods without 
jumps. All jumps above half the local standard deviation of the 
signal-removed TOI are found. 

3.9. Flagging samples 

We now build a total flag out of several flags in order to qualify 
the TOI for mapmaking. A sample with any of the following 
flags is considered invalid data: 

1. The unstable pointing flag described in Sect. 3.1 (see also 
Planck Collaboration I 2013), typically accounting for three 
minutes per pointing period (roughly seven percent of data). 

2. The missing or compression error data flag, discard- 
ing a fraction of 10"^ of the whole mission (see 
Planck HFI Core Team 201 la). 

3. The bolometer plate temperature fluctuation flag (see 
Sect. 3.5). One to two percent of the data are flagged this 
way depending on time. 

4. The glitch flag. Typically between eight and twenty percent 
of data are flagged depending on the time and the bolometer 
(Planck Collaboration X 2013). For PSBs, both detectors are 
flagged even if, in a few cases, only one of them sees the 
glitch. 

5. The jump flag. Two hundred samples per jump are flagged, a 
fraction less than 10"^ (see Sect. 3.8). 

6. The Solar System Object (SSO) flag. A zone of exclusion 
is defined around Mars, Jupiter, Saturn, Uranus, Neptune, 
and around the 24 asteroids detected at 857 GHz (see 
the complete discussion of asteroids in Appendix A of 
Planck CoUaboration XIV (2013)). The coordinates of the 
objects are obtained from the JPL Horizons database^ 
(Chamberlin et al. 1997), which uses the actual position of 
the Planck Satellite in its orbit around the L2 point. As the 
Solar System Objects are moving in celestial coordinates, 
the data (although valid) must be discarded and this masking 
must be done in the time domain. One survey map has thus 
up to 35 holes, which are filled by information from other 
surveys. Overall, the final maps have almost no holes: less 
than 10"^ of pixels are missing. 

On top of flagged samples in the TOI, we show now that some 
flagging also has to be done at the pointing period level (i.e., 
considering entire rings of data). 

3.10. TOI Qualification 

We will later make maps by projecting the PBR into sky co- 
ordinates. We therefore base our criteria for the acceptance or 
rejection of a given pointing period for mapmaking by using the 
statistics of only the signal-removed TOI as an estimate of the 
noise. Hence, the qualification is minimally biased by the signal 
itself. 

This noise estimate timeline is first analysed in the Fourier 
domain. In Planck HFI Core Team (2011b) we described how 
the detector white noise level (equivalently, the NET) is de- 
termined from noise periodograms in the frequency region bey- 
ond the low frequency excess noise component and before the 
increase introduced by the time-constant deconvolution. This 
white noise level differs from that deduced directly from the 
maps. We have therefore adopted a method for determining the 

^ http://ssd.jpl.nasa.gov/ 
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"total" noise: this is derived directly from the standard deviation 
of the signal-removed TOI after flagging the Galactic plane and 
point sources for each pointing period. The standard deviation 
is then corrected upward by a term that depends on the duration 
of the pointing period. This term has the form a/7^7(7^^^' 
where d is the pointing period duration in minutes and / is the 
fraction of valid data within a pointing period (typically 0.8). 
It accounts for the fact that a fraction of the noise remains in 
the PBR, and that fraction is smaller on longer pointing peri- 
ods, where the signal is better estimated. The final value of this 
total noise level for each detector (third column of Table 1) is 
then taken to be the peak of the distribution of the root-mean- 
square (rms) noise of the valid pointing periods. The typical bo- 
lometer NETs, measured on the deconvolved TOI between 0.6 
and 2.5 Hz, are also given in Table 1 (second column). The total 
noise can be viewed as proportional to the square root of the 
integral of the noise-equivalent temperature across the total fre- 
quency bandpass (including filtering eff'ects) from to 91 Hz. It 
is normalised to one second of integration time so that pure white 
noise would have numerically identical NET and total noise. 
Deconvolution and filtering eff'ects can produce a NET larger or 
smaller than the total noise. 



Table 1. Noise characteristics. The quoted value is the noise that 
is obtained on the map after one second of integration, i.e., about 
180 hits for one bolometer. This gives the noise characteristics 
of bolometers, inverse quadratically averaged over the similar 
ones within a channel ("P" is for polarization sensitive bolomet- 
ers and "S" for spider web unpolarized bolometers). The second 
column gives the white noise level measured on the power spec- 
trum in the range of 0.6 to 2.5 Hz. The third column gives the 
total noise, i.e., the rms noise at the map level for one second 
of integration time of a single detector. The next column recalls 
the goal stated before launch by Lamarre et al. (2010). Note that 
the mean integration time per detector per (1.7')^ pixel is 0.56 
seconds (0.63) for the 100-353 (resp. 545-857 GHz) channels, 
for the nominal mission (see Fig. 21). 



Band 


NET 


Total Noise 


Goal 


Units 


lOOP 


71 


132 


100 


A/KcMB S^/^ 


143P 


58 


65 


82 


^KcMB S^^^ 


143S 


45 


49 


62 


yUKcMB S^/^ 


217P 


88 


101 


132 


//KcMB S^^^ 


217S 


74 


66 


91 


ywKcMB s^/^ 


353P 


353 


397 


404 


^KcMB S^^^ 


353S 


234 


205 


277 


A^KcMB S^^^ 


545S 


0.087 


0.052 


0.116 


MJy sr-1 s^/2 


857S 


0.085 


0.056 


0.204 


MJy sr-i s^/2 



The noise estimate is also analysed in the time domain. 
Figure 4 shows one pointing period of the signal-removed TOI 
for four detectors along with the histogram of the samples. 

Figure 5 shows an example of the evolution of this total noise 
(uncalibrated at this stage) during the first four surveys (hence 
beyond the nominal mission). The bias linked to the pointing 
period duration has been corrected for. Notice that the total 
noise is mostly constant with peak-to-peak variations at the per- 
cent level, suspected to be induced by ADC non-linearities (see 
Sect. 6.2). 



A pointing period that deviates significantly from the others 
is flagged in the qualification process and not used in further 
processing. Three criteria were used to quantify the deviation: 

- the absolute value of the diff'erence between the mean and 
the median value of the signal-removed TOI, 

- the total noise value, and 

- the Kolmogorov-Smimov test deviation value. When a given 
pointing period is stamped as anomalous for more than half 
the bolometers, it is discarded for all of them. 

The common anomalous pointing periods are linked to some 
spacecraft events or strong solar flares, while the individual an- 
omalous pointing periods correspond to incidents in the overall 
level of the timeline (like strong drifts). 

There are two detectors (one each at 143 and 545 GHz, not 
counting toward the total of 50 HFI detectors) that are not used at 
all because they show a permanent non-Gaussian noise structure, 
which we denote random telegraphic signal (RTS). The RTS bo- 
lometer signal-removed TOI show an accumulation of values at 
two to five discrete levels. The measured signal-removed TOI 
jump at random intervals of several seconds between the dif- 
ferent levels. The jumps are uncorrected. A third detector (at 
857 GHz) shows the RTS phenomenon but only for some well- 
delimited periods. Other detectors show a significant but smaller 
RTS in episodes of short duration. A dedicated module searches 
for these episodes, which are then added to the list of discarded 
pointing periods. Ten pointing periods of exceptionally long dur- 
ation are discarded for all bolometers, because noise stationarity 
is not satisfied in those cases. 

The discarded pointing periods represent less than 0.8 % of 
the total integration time for the nominal mission. 

Many examples of the characterization and qualification of 
the TOI processing can be found in Planck Collaboration ES 
(2013), where simulations are used to limit the possible eff'ects 
of undetected RTS periods to a fraction of the overall noise level. 

In addition, the explanatory supplement presents a number 
of simulations and other tests which have been done to set lim- 
its on possible contamination from a number of other possible 
systematic effects. 

3.11. Big planet TOI 

As the TOI of the standard pipeline are interpolated around 
Mars, Jupiter, and Saturn, they cannot be used for the beam re- 
construction (see Sect. 5). Therefore a dedicated pipeline is run 
on appropriate pointing period ranges. The diff'erence with the 
standard pipeline is mainly a special deglitching step done on 
the planet samples to cope with such a strong signal and its vari- 
ations during a given pointing period. Also, the jump correction 
is not applied to avoid its triggering by the planet transit. 

4. Detector pointing 

Here we summarize the pointing solution we use, determined 
with the overall goal of limiting detector pointing errors to 
less than 0.'5 (rms). The satellite pointing comes from the star 
tracker camera subsystem, which gives the location of a fidu- 
cial boresight direction as a function of time, sampled at 8 Hz. 
In practice, science data analysis requrires the pointing of each 
detector, sampled at the same rate as the detector signal. This, in 
turn, requires a number of separate steps: 

1. Resampling to the acquisition frequency (180.4 Hz); 
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Figure 4. The signal-removed TOI sample values are displayed for one pointing period and four detectors (left panel of the four 
sets). Flagged samples are not shown. Samples falling near the Galactic plane or point sources are not shown either (corresponding 
to one-minute periodic blanks in the plots). The right panel of the windows shows the histogram (in black) of the samples. The 
coloured curve is the symmetric version around the mean. 



2. Rotation from the fiducial boresight to the detector line of 
sight; and 

3. Correction for detector aberration 
(Planck Collaboration XXVII 2013). 

The second step above must be calibrated in situ: the detailed 
geometry must be measured by comparison with the positions 
of known objects (planets and other bright sources whose ex- 
tent is small compared to the width of the beam). The primary 
source of this geometrical calibration is the planet Mars, which 
is bright (but not so bright as to drive the detectors into a non- 
linear response) and nearly pointlike (with a mean disk radius of 
4!'l 1 1 during the nominal mission). Figure 6 shows the pointing 
of each detector relative to the pre-launch optical model. 

Note that the in- scan pointing is degenerate with any phase 
shift induced by the combined detector time response and op- 
tical beam and any attempt to remove them by deconvolution. 
For that reason, this calibration of the pointing cannot itself be 
considered as a measurement of the focal plane geometry per se. 

Because the transfer of the pointing from the star tracker to 
the detectors depends on the satellite rotation axis and hence the 
moment of inertia tensor, we do not expect the detector pointing 
relative to the star tracker to be constant in time. We therefore 
reconstruct the pointing solution relative to the star tracker as a 



function of time, and in the following discuss the accuracy of the 
reconstruction in diff'erent frequency regimes. 

High-frequency fluctuations in the spacecraft and focal plane 
orientation are inhibited by inertia but are still present in the raw 
star tracker pointing. Having no useful signal content, this re- 
gime is treated by applying a low pass filter to the reconstruc- 
ted pointing. Similarly, high-frequency errors can be assessed 
by studying the noise-dominated high frequency pointing power 
spectrum. 

Intermediate-frequency pointing errors, on one minute time 
scales, are characterized with Jupiter transits. We make use of 
the optimal combination of high signal-to-noise ratio and rel- 
atively wide beam at the Planck 143 GHz channel and first do 
a global fit of the planet position in the transit data. We then 
consider a small pointing off'set to each 60 second scanning 
circle. The test succesfully recovers an expected interference sig- 
nal from the radiometer electronics box assembly (REBA; see 
Planck Collaboration II 2011 and Planck Collaboration I 2013). 
Thermal control on REBA was adjusted on day 540 after launch 
and subsequently Jupiter transit three does not exhibit the inter- 
ference. We find that the intermediate-frequency pointing error 
is less than y (rms) before the REBA adjustment and less than 
V beyond day 540 (see Fig. 7). 
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Figures. Top: Total noise estimated for each stable pointing period for a 143 GHz PSB. A smoothed version is shown in green . 
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Figure 6. Position of the individual HFI detectors with respect 
to a pre-launch model as measured on the first passing of Mars. 
The horizontal is the direction of the scan ("in-scan"); vertical is 
perpendicular to this ("cross-scan"). 



The pointing is easiest to monitor over long time scales, us- 
ing observations of bright planets. In Fig. 8 we show the dif- 
ference between the first two observations of Mars; diff'erences 
along the scan direction have a mean of 0''85 and a standard de- 
viation of 0759', cross-scan diff'erences in the direction perpen- 



dicular to this are 371 ± 27^ and are systematic with frequency. 
Indeed, analysis of further planet-crossings and of aggregate 
high-frequency point-source crossings exhibits arcminute- scale 
evolution over the course of the mission. Fitting for an overall 
off'set for each planetary transit allows us to measure long time 
scale (several days) pointing fluctuations. We complement this 
sparse sampling of the overall ofl'set by time domain position 
fits of brightest compact radio sources as well as monitoring of 
many low-flux high-frequency sources in aggregate. In the pro- 
cess we learned that each radio source position was biased by a 
few arcseconds from their catalogued position. We take the bias 
to result from convolving the source emission with our effect- 
ive beams and correct for it by assuming that the low frequency 
pointing correction derived from planet transits alone is already 
accurate enough for debiasing the positions. For each compact 
source we solve for its apparent position in the Planck data by 
considering a single constant ofl'set to its catalogued position. 

The planet positions alone indicate a change in the average 
level of the low frequency fluctuations between Jupiter transit 
two and Neptune transit three (days 418 and 540). The ra- 
dio sources, particularly QSO J0403-3605, further narrow the 
transition before day 456. It is likely that the ofl'set is due 
to thermoelastic deformation caused by the sorption cooler 
system switchover on day 460 (Planck Collaboration II 2011; 
Planck Collaboration I 2013). If entirely untreated, the low fre- 
quency pointing error dominates the error budget with around 
15'' (rms) error. Two independent focal plane solutions bring the 
low-frequency rms error well below 10'' and applying a time- 
dependent correction reduces the low frequency error to the or- 
der of few arcseconds. 

Figure 9 shows the uncorrected planet positions with the 
debiased bright point- source positions along with the resulting 
Planck pointing solution, which smoothly interpolates with a 
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Figure 7. Measured intermediate frequency pointing errors at one minute intervals. The offsets were smoothed with a five-point me- 
dian filter. Vertical lines mark the pointing period boundaries. Overplotted solid line is the scaled and translated REBA temperature. 
REBA thermal control was adjusted between Jupiter transits two and three. Left: Pointing fluctuations during the second Jupiter 
transit. Right: Fluctuations during the third Jupiter transit, after the REBA thermal control was adjusted. Notice the diff'erent vertical 
scales. 



spline fit between planet observations. Figure 8 shows the off'- 
set of Mars for all detectors between the first and second surveys 
after correction by the pointing solution and Fig. 10 the eff'ect 
of the corrected pointing on the high-frequency point sources at 
545 and 857 GHz. 

Afer these corrections, the pointing error is about T!5 (rms) 
over the mission, varying from less than V during planet ob- 
servations, and increasing between planet crossings. The indi- 
vidual HFI detector locations (along with their respective beam 
patterns) are show in Fig. 11. 
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Figure 8. The diff'erence between the pointing for the first and 
second crossings of Mars, as in Figure 6, but note the diff'erent 
axis scales. 



5. Detector beams 

Here we provide a brief overview of the measurement of the 
Planck beams and the resulting harmonic- space window func- 
tions that describe the net optical and electronic response to the 
sky signal (for a full discussion see Planck Collaboration VII 
2013). 

At any given time, the response to a point source is given 
by the combination of the optical response of the Planck tele- 
scope (the optical beam) and the electronic transfer function. 
The latter is partially removed as discussed in Sect. 3 and 
Planck Collaboration VII (2013) but any residual is taken into 
account as part of the beam. This response pattern is referred to 
as the scanning beam. However, the mapmaking procedure dis- 
cussed in Sect. 6 below means that any map pixel is the sum of 
many diff'erent elements of the timeline, each of which has hit 
the pixel in a diff'erent location and in a diff'erent direction. Thus, 
the effective beam (Mitra et al. 2011) takes into account the de- 
tails of the scan pattern (and even this is a somewhat simplified 
picture, as long-term noise correlations mean that even samples 
taken when the telescope was looking elsewhere contribute to 
a given pixel). Finally, the multiplicative eff'ect on the angular 
power spectrum is encoded in the effective beam window func- 
tion (Hivon et al. 2002), which include the appropriate weights 
for analyzing aggregate maps across detector sets or frequencies. 

5.1. Scanning beams 

For the single-mode HFI channels, the scanning beam is 
well-described by a two-dimensional elliptical Gaussian with 
small perturbations at a level of a few percent of the peak 
(Huff'enberger et al. 2010). We measure the beams on obser- 
vations of Mars, and model them with B- Spline fits to the 
timeline supplemented with a smoothing criterion as discussed 
in Planck CoUaboration VII (2013). 

We measure the beam on a square patch 40' on a side. This 
main beam pattern must be augmented by several eff'ects that are 
visible within roughly one degree of the beam centre due to vari- 
ations of the mirror surface. At all frequencies we see evidence 
of un-modelled eff'ectively random dimpling, well-fit by a sum 
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Figure 9. Compilation of planet and bright point-source positions versus time since the start of the mission. We denote the REBA 
thermal control adjustment on day 540 with a dashed grey line and the end of the nominal mission (day 563) with a solid grey line. 
Note the two discontinuities added to the pointing correction model at days 455 and 818. The dashed green line denotes the spline 
fit to the planet positions used to correct the long-term pointing variations. Left: Compilation of measured planet position off'sets 
across the Planck frequencies. Usable planets change by frequency: the blue spectral energy distribution of the planets renders all 
but Jupiter too dim for required positioning in the LFI frequencies and Jupiter and Saturn positions above 217 GHz are compromised 
due to nonlinear bolometer response. Right: Debiased bright source position off'sets in Planck lOOGHz-217 GHz data compared 
to a trend line fitted to the planet position off'sets. Each source was fitted for an apparent position rather than using the catalogued 
position to measure pointing off'set. Two successive observations of the same source have opposite scanning directions so having 
the wrong apparent position produces opposite pointing errors that we have corrected for. 



of Ruze (1966) models. At 353 GHz and higher we see evidence 
of the hexagonal backing structure of the mirrors, included dir- 
ectly in the B- Spline fits. The scanning beam patterns for each 
detector are shown in Fig. 11. 

Other steps in the data processing pipeline also aff'ect the 
scanning beam pattern. We see residual eff'ects from the decon- 
volution as a shoulder in the beam localised to the "trailing" side 
of the scan. The residual pointing uncertainty results in a slight 
broadening of the beam which is modelled in the Monte Carlo 
simulations we use to estimate the beam errors. Further simula- 
tions indicate that undetected or inaccurately removed glitches 
(Sect. 3.3) result in a negligible additional bias, as do the uncor- 
rected ADC non-linearities. Because the response of the detector 
depends on the spectral shape of the source, there is a colour 
correction (see Sect. 6.6) from the approximately known planet 
spectrum (roughly y^) to the nominal y"^ shape. 



All of these give a change in the beam solid angle of roughly 
0.5 % at 100 GHz and less than 0.2 % for 143 GHz and above. 
However, there is an additional contribution to the scanning 
beam from near sidelobe response beyond the square patch on 
which the beam is measured (e.g., from the deconvolution of the 
time response). Simulations show that the near sidelobes con- 
tribute less than 0.2 % of the beam power. 

5.2. Effective beams 

These scanning beams are used to calculate the eff'ective beam 
response at a given pixel. A symmetric beam in a uniformly- 
sampled pixel (and ignoring eff'ects induced by long-time- scale 
noise correlations accounted for in the mapmaking procedure) 
would give an eff'ective beam that is the convolution of the 
scanning beam pattern with a top-hat pixel window function. 
Departures from these idealities due to scanning beam asym- 
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Figure 10. Pointing evolution from aggregate high-frequency point sources, measured by comparing point sources seen at 545 and 
857 GHz to known IRAS positions. Light-colored points show individual source deviations, points with error bars give ten-day 
errors. Different colors correspond to the two frequencies and four individual Planck sky surveys. 



metry and the scan strategy means that the effective beam will 
be asymmetric and depend upon sky location: ignoring long- 
term noise correlations, the effective beam is given by averaging 
the scanning beam pattern at the observed locations and orienta- 
tions of the real scan pattern in each pixel. In practice, this must 
be approximated by using a coarse pixelization of the scanning 
beam or restricting the calculation to those components with low 
spherical-harmonic multipoles m. We have developed a set of 
methods in real and harmonic space to account for these effects, 
described more completely in Planck Collaboration VII (2013) 
and references therein. 

We propagate the scanning beam pattern from each detector, 
along with the largest- variance eigenmodes of the Monte Carlo 
error covariance matrix, using the Planck scanning strategy. The 
error eigenmodes are scaled by a factor of 2.7 to account for the 
unmodelled near sidelobe bias in the scanning beam. The even- 
tual calculation of power spectra requires the cross-correlation 
of various detector pairs and more generally arbitrary weighted 
combinations of detector pairs used to construct the final CMB 
power spectrum. Our algorithm results in both mean beam pat- 
terns (Fig. 12) and sky-averaged window functions, along with 



error eigenmodes on the latter (Fig. 13), for the required detector 
combinations. The far sidelobe response of the instrument (re- 
sponse to signal more than five degrees from the beam centroid, 
discussed in Sect. 6.4) negligibly biases the calibration and 
the effective beam window function (Planck Collaboration VII 
2013; Planck Collaboration XIV 2013). The full window func- 
tions are available with the Planck data release, and the key nu- 
merical results are shown in Table 4. 



6. Mapmaking and photometric calibration 

This section gives an overview of the mapmaking and pho- 
tometric calibration procedure. More details are described in 
Planck Collaboration VIII (2013). 

6.1. Map projection & calibration techniques 

The HFI mapmaking scheme is very similar to that described 
in Planck HFI Core Team (201 lb). The inputs for each detector 
come from the TOI processing (see Sect. 3) with the associ- 
ated invalid data flags. Mapmaking is done in three steps. We 
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Figure 11. The scanning beam patterns reconstructed from planet observations at their respective positions in the HFI focal plane. 
The beams are plotted in logarithmic contours of -3, -10, -20 and -30 dB from the peak. Polarization sensitive bolometers are 
indicated as a pair of contours, one black and one blue (in most cases the underlying black contour is invisible, indicating essentially 
identical beam shapes). 



first take advantage of the redundancies during a stable point- 
ing period (ring of data) by averaging each detector's measure- 
ments in HEALPix (Gorski et al. 2005) pixels, building a new 
data structure, the HEALPix pixel ring (HPR). We prefer to use 
this structure for mapmaking, rather than the phased binned rings 
previously described, because the latter introduce an additional 
smoothing in the reconstructed map power spectra. We set the 
HEALPix resolution for the HPR at the same level as for the 



final maps, A^side = 2048, for the same reason. In the second 
step we use these HPR to perform the photometric calibration of 
each detector's data. Once this is done, the mapmaking per se 
is performed. As for the previous release, the mapmaking pro- 
cedure algorithm we used is an implementation of a destriping 
algorithm, which is made possible by the characteristics of the 
HFI detector noise. 
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Figure 12. The effective beam solid angle (upper panel) and the 
best-fit Gaussian ellipticity (lower panel) of the 100 GHz eff'ect- 
ive beam across the sky in Galactic coordinates. 



In the destriping framework, detector noise is modelled as a 
set of constants, called off'sets, representing the low-frequency 
drift of the signal baseline over given time intervals, and a white 
noise component uncorrelated with these constants. The general 
mapmaking scheme may thus be written in this approximation 
as 



d = GxAT + ro + n, 



(2) 



where T represents the pixelized sky [which may be a vector 
(/, 2, U) Sit each pixel if polarization is accounted for], G is the 
detector gain, A is the pointing matrix indexed by sample num- 
ber and pixel, equal to zero when the pixel is unobserved at that 
sample, one when the pixel is observed by an unpolarized de- 
tector, or the vector (1, cos 2a, sin2Qf) for a polarized detector 
at an angle a with respect to the axis defining the Stokes para- 
meters, and r is the matrix folding the ring onto the sky pixels. 
From the above equation, the off'sets o are derived through max- 
imum likelihood, imposing an additional constraint, in our case 
that the sum of the off'sets has to be equal to zero (arbitrarily, 
since we do not measure the absolute temperature on the sky, 
only diff'erences). The performance of this implementation has 
been evaluated using simulations in Tristram et al. (201 1). 

We produce temperature (and polarization whenever pos- 
sible) maps for a number of data sets: 

- individual detectors, 

- independent detectors sets at the same frequency (see 
Table 3), 

- the combination of all detectors at each frequency. 
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Figure 13. Eff'ective beam window functions (solid lines) for 
each HFI frequency. The shaded region shows the full ±lcr er- 
ror envelope. Dashed lines show the eff'ective beam window 
function for Gaussian beams with full-width at half-maximum 
(FWHM) parameter 9.'65,7.'25,4.'99,4.'82,4.'68, and 432 for 
100, 143, 217, 353, 545, and 857 GHz respectively. We com- 
pute the 100 GHz window function to ^ = 2500, 143 GHz to 
^ = 3000, and the rest to ^ = 4000. 



To enable systematic checks, we build maps for the whole time 
interval spanned by the HFI data as well as more restricted in- 
tervals corresponding to individual surveys. We also build maps 
in each of the above cases, using data from the first and second 
half of the rings independently. In each case (a full ring or its 
two halves) we first determine the off'sets on the total time range 
and apply them to produce maps on any more restricted time 
range. Together with each temperature map produced, we build 
other maps containing respectively the hit count (number of TOI 
samples in each pixel) and the propagation of the TOI sample 
variance, as measured by the total noise NET (see Sect. 3.10). In 
the end, a total of more than 6500 maps are produced. 

6.2. Absolute photometric calibration 

We use two techniques to perform the photometric calibration 
of the HFI detectors, one for lower frequency channels (100- 
353 GHz, typically given in temperature units, Kcmb), another 
for higher frequency channels (545 and 857 GHz, given in flux 
density units, MJy sr"^). In both cases, significant changes have 
occurred with respect to the early data processing. The photo- 
metric calibration processes and their performance are described 
in detail in Planck Collaboration VIII (2013). Our philosophy is 
to calibrate each detector individually, without relative calibra- 
tion between them. 

Previously (for the early results from Planck), the lower fre- 
quencies were calibrated using the Solar dipole, and the higher 
frequencies using FIR AS dust spectra measurements assuming 
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constant gain. When more data were accumulated, comparisons 
between measurements taken one year apart in identical detector 
configuration unambiguously showed that the detectors' gains 
show variations of the order of 1 to 2 %, over a wide range of 
time intervals (one day to one year). 

We showed that these apparent gain variations are first order 
consequences of non-linearities of the ADC used in the read- 
out electronic chain (see Sect. 7.2.1). Correcting for such eff'ects 
requires a precise knowledge of the ADC transfer function. To 
acquire these we characterize the read-out response using warm 
data since the end of the 100 mK cooling in January 2012. This 
data- taking was not completed in time for the 2013 HFI data 
release. At the frequencies where the Solar dipole has a high 
enough amplitude (100 to 217 GHz) we evaluated the apparent 
gain variations of the detectors by solving the non-linear equa- 
tion data = gain x sky -h noise for both gain and sky signal. 
Examples of gain- variation measurements are show in Fig. 14. 
Note the wide range of behaviour, from slow drifts (e.g., 100- la) 
to apparent jumps (e.g., 143 -la). 
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Figure 14. Samples of relative gain variations reconstructed for 
six CMB channel bolometers. Data for each bolometer have 
been translated vertically for clarity, with a 3 % spacing. Thin 
lines show the averaged levels in each case. The observed vari- 
ations range in general over ±1.5 %. 



tematics) after correction are lower than 0.3 %, as measured by 
comparison with the Solar dipole. 

The presence of residual non-linear systematic biases in our 
data precludes the use of potentially more precise techniques, 
such as those discussed in Tristram et al. (201 1) using the orbital 
CMB dipole anisotropy. In Planck Collaboration VIII (2013) we 
show that using a calibration derived from the orbital dipole 
would lead to larger detector-to-detector relative calibration dis- 
persion, which induces large-angular-scale patterns in polariza- 
tion through dipole leakage. In addition, for detectors for which 
the Solar and orbital dipole diff'erence is large (typically around 
0.5 % but in few cases as large as 1 %), a residual dipole is appar- 
ent after subtracting the WMAP prediction. Thus, the calibration 
for 100 to 353 GHz is based on the Solar dipole as measured by 
WMAP (we use a non-relativistic calculation; the roughly 0.1 % 
relativistic correction is smaller than other residuals such as the 
uncorrected ADC non-linearities). 

For the two highest frequency channels, the calibration 
scheme that was used for the early data release was based on 
FIRAS data. As described in Planck Collaboration VIII (2013), 
several studies Since the first release indicated that the calibra- 
tion scale of these two channels was somewhat incorrect. These 
include measurements of the dust spectral energy distribution 
and of the amplitude of the CMB anisotropy itself (still detect- 
able at 545 GHz). An unpublished re-analysis of the FIRAS data 
(Liang et al. 2012) also seems to indicate that the FIRAS spectra 
might contain unexpected additional systematic errors. Since we 
also observed systematic diff'erences between HFI and FIRAS 
data, we decided to revise our calibration strategy. We therefore 
now supplement our calibration using the reconstructed fluxes of 
planets, compared with models of their emission (Moreno 2010), 
in order to renormalize the relative FIRAS calibration at 545 and 
857 GHz. We used observations of Uranus, Neptune, and Mars 
(Jupiter and Saturn, although observed in Planck data, are too 
bright and hence afl'ected by detector non-linearities). Compared 
to previous data releases, this new calibration scheme amounts to 
a decrease in flux densities of 15 % and 7 % at 545 and 857 GHz, 
respectively. 

Various techniques have been used to evaluate the relat- 
ive (intra- or inter-frequency) and absolute calibration accur- 
acy (Planck Collaboration VIII 2013). Results are summarized 
in Table 4. The main limitations on the photometric calibration 
accuracy that we have identified result from the residual non- 
linearities for low frequency channels, and from the accuracy of 
the models of planetary brightness at 545 and 857 GHz. 

For this data release, the zero levels of the maps, which 
Planck cannot determine internally, have not been set. Still, 
we estimate in (Planck Collaboration VIII 2013) a Galactic and 
an extragalactic zero level (Tables 4 and 5). For the Galactic 
case, we compute the map brightness corresponding to zero gas 
column density as traced by the 21 cm emission from neutral 
hydrogen. For the monopole term of the CIB, we use an empir- 
ical model, which can be used to enable total emission analysis. 
These estimates can be used to set the appropriate zero levels of 
HFI maps and are also summarised in Table 4. 



These measurements have been used to correct the bolo- 
meter data in the present data release. As these gain variations 
are just the first-order consequence of the ADC non-linearities, 
these corrections will not remove all ADC non-linearity eff'ects. 
We estimate (Planck Collaboration VIII 2013) that the remain- 
ing apparent gain variations (or other residual time- variable sys- 



6.3. Overview of HFI map properties 

Figures 15 to 20 show the sky maps constructed from HFI data. 
These correspond to maps uncorrected for zodiacal light emis- 
sion and far side lobe pick-up (discussed in the next section). 
The top row of the figure series shows the intensity maps re- 
constructed from the nominal mission, while the second row 
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shows the difference between maps made of the first and the 
second half of each stable pointing period (i.e., half-ring maps), 
which provides intuitive evidence of the level and distribution of 
the residuals, since these maps illustrate the diff'erence between 
maps constructed from rings taken about 20 minutes apart; there- 
fore they illustrate the variation of the detector signal along a 
sky circle on such a time-scale, which certainly dominates the 
rms. Longer time- scale variations can be judged from the third 
row of plots showing the diff'erence in the sky between the first 
and second surveys. This time one sees variations on a six- 
month time scale, which maximizes systematic eff'ects by com- 
paring measurements taken with a diff'erent scan orientation of 
the satelitte. Note that for 100-217 GHz, the colour scale of the 
second row is enlarged by a factor of 75, by 200 at 353 and 
545 GHz, and by a factor greater than 730 at 857 GHz. 

The distribution of integration time per channel is shown in 
Fig. 21. The map noise properties may be evaluated from vari- 
ance maps computed from the integration time per pixel, but a 
better approach uses the half-ring map diff'erences (middle row 
of Figs. 15 to 20), which gives a better rendition of the noise at 
map level due to the non-white structure of the timeline noise 
after time-response deconvolution. We show the power spectra 
of these diff'erences, for each frequency, in Fig. 22. With respect 
to the Planck early data release, these noise spectra are notably 
flatter at high £ thanks to the use of a better low pass filter in the 
TOI processing. The figure also shows the spectra of the diff'er- 
ences between maps made from the first and second sky surveys 
(bottom row of Figs. 15 to 20), giving an indication of the con- 
tribution to the noise from longer time scales. The two sets of 
spectra converge at high L 

The average of the half-ring power spectra from f = 100 to 
6000 are reported in Table 4. 

One should note, however, that the half-ring maps give an 
estimate of the noise that is biased low (by a couple of per- 
cent) due to small correlations induced by the way the timelines 
have been degliched. This is discussed in detail with the help of 
simulations in Sect. 7.2.2. Estimates of the noise level from the 
timelines and from the maps were discussed in Sect. 3.10 and 
given in Table 1. They are consistent with alternative descrip- 
tions of the noise which were computed in the HFI calibration 
and map-making paper (Planck Collaboration VIII 2013). The 
summary Table 4 (lines cl-c2) gives sensitivity levels derived 
from these analyses. 

6.4. Far-side lobe and zodiacal light correction 

Most of the signals Planck observes on the sky are fixed to the 
celestial sphere. There are, however, two exceptions: zodiacal 
light emission (ZLE) and far side lobe (FSL) contamination. 
Because Planck surveys the sky twice per year with a "cycloidal" 
scanning strategy (Dupac & Tauber 2005), each time a given 
location on the (distant) sky is observed, a diff'erent column of 
nearby dust in our Solar System is sampled, leading to a slightly 
diff'erent ZLE. The signature of zodiacal emission is indicated 
by the "S"-shaped band along the Ecliptic plane seen in the 
857 GHz diff'erence between maps made six months apart (bot- 
tom row of Fig. 20), when Planck is located at antipodal points 
in its orbit and therefore looking through very diff'erent columns 
of zodiacal dust. Figure 23 shows the sky pattern of the ZLE 
model in each of the HFI bands. 

In addition, there is some small, asymmetric, oflf-axis beam 
response that also leads to additional time variability. This is 
seen as the arcs near the north and south Galactic poles in 



Fig. 20. Both of these are most pronounced in the highest- 
frequency channels. We remove these signals so that further 
Galactic and extragalactic studies will not suff'er from interplan- 
etary dust emission. 

The removal process is described in more depth in 
Planck Collaboration XIV (2013), and summarized here. It con- 
sists of fitting survey maps from each horn to the COBE zodi- 
acal emission model (Kelsall et al. 1998) to find the emissivities 
at the HFI wavelengths, plus templates for the Galactic and di- 
pole signal seen through a model of the far sidelobes. With the 
amplitudes for each template from this fit, we again turn to the 
COBE and far sidelobe models to reconstruct the implied zodi- 
acal emission at all times during the survey for each individual 
detector. We convert the values to yuKcMB, and remove these sig- 
nals before the data are recombined into maps. 

The fit coefficients for the zodiacal emission are presented 
in Planck Collaboration XIV (2013), as are representative maps 
of the amount of signal removed from each map, and the power 
spectra of these maps. At 143 GHz, the zodiacal emission con- 
sists essentially of a band of emission between Ecliptic latitudes 
of roughly ±15°, reaching a maximum of 5yuKcMB- The power 
spectrum at 143 GHz is small: it appears only at multipoles less 
than roughly 75, and expressed as (InY^ ^(^ -h 1)Q is every- 
where less than 2.5 /^K^^g in all multipoles (and decreased by 
binning). 

Maps are made available both with and without ZLE re- 
moval. In addition, as there is a detection of the Galaxy seen 
through the far sidelobes in some horns at 857 and 545 GHz, 
these signals are removed as well, for these frequencies only, 
along with ZLE removal. Figure 23 shows the sky pattern correc- 
ted for, and Fig. 24 shows its power spectrum for 100-353 GHz, 
separately for odd and even multipoles, reflecting the approx- 
imate North/South symmetry of the pattern. The correction is 
much less uncertain at the three highest frequencies, but the 
low frequency spectra allows to gauge the rough level of the 
contribution to the low-^ power spectra in uncorrected maps. It 
shows that in the moderate-^ part of the Planck likelihood, i.e., 
at ^ > 50, the ZLE contribution in uncorrected maps is smaller 
than 1 /^K^MB and can thus be neglected. But in order to extract 
the CMB, we found that one cannot use only a partial correction 
of the ZLE at the frequencies where it is well determined, i.e., 
at 353 GHz and higher frequencies, since lower- amplitude "S"- 
shaped residuals (in the CMB map) are found when the correc- 
ted maps are used. Nevertheless, overall, we discovered that the 
best CMB cleaning was achieved by the full multi-component 
CMB and foreground separation process itself on maps not oth- 
erwise explicitly corrected for the ZLE. This is addressed in 
more detail in Planck Collaboration XII (2013), but the main 
reason can be seen right away: the relative weights of the maps 
that are used to generate a CMB map are found to be the same 
when one set of HFI maps or the other is used, and are given 
by (30GHz:-0.01, 40 GHz: -0.04, 70 GHz: -0.16, 100 GHz: 
-0.22, 143 GHz: 1.2, 217 GHz: 0.33, 353 GHz: -0.11, 545 GHz: 
3x10""^, 857 GHz: 1.8x10"^ for the NILC component- separation 
technique described in Planck Collaboration XII 2013). The typ- 
ical peak-to-peak amplitude of the zodiacal correction is about 
(70 GHz: 6, 100 GHz: 6, 143 GHz: 5, 217 GHz: 15, 353 GHz: 
150, 545 GHz: 2400, 857 GHz: 2 x one sees that with 

the above weights the 353-857 GHz contribution roughly can- 
cels the 143-217 GHz one, while lower frequency channels have 
both a lower zodiacal contribution and a lower weight in determ- 
ining the CMB map. 
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Figure 15. HFI maps at 100 GHz. The first row gives the intensity in yuKcMB- The second row shows the diff'erence between maps 
made of the first and the second half of each stable pointing period (i.e., half-ring maps). The third row shows the diff'erence between 
the first and second survey maps. 
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Figure 16. HFI maps at 143 GHz. The first row gives the intensity in yuKcMB- The second row shows the diff'erence between maps 
made of the first and the second half of each stable pointing period (i.e., half-ring maps). The third row shows the diff'erence between 
the first and second survey maps. 
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Figure 17. HFI maps at 217 GHz. The first row gives the intensity in yuKcMB- The second row shows the diff'erence between maps 
made of the first and the second half of each stable pointing period (i.e., half-ring maps). The third row shows the diff'erence between 
the first and second survey maps. 
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Figure 18. HFI maps at 353 GHz. The first row gives the intensity in yuKcMB- The second row shows the diff'erence between maps 
made of the first and the second half of each stable pointing period (i.e., half-ring maps). The third row shows the diff'erence between 
the first and second survey maps. 
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Figure 19. HFI maps at 545 GHz. The first row gives the intensity in MJy sr"^ The second row shows the diff'erence between maps 
made of the first and the second half of each stable pointing period (i.e., half-ring maps). The third row shows the diff'erence between 
the first and second survey maps. 
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Figure 20. HFI maps at 857 GHz. The first row gives the intensity in MJy sr" . The second row shows the diff'erence between maps 
made of the first and the second half of each stable pointing period (i.e., half-ring maps). The third row shows the diff'erence between 
the first and second survey maps. 
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Figure 21. HFI channel integration time distribution function. 
For each HFI channel, the histogram of the integration time t 
per pixel of area (1.7')^ is shown. The black curve is for the 
nominal all-sky (/sky = 1) mission. The red curve represents 
the masked sky, using the CL39 (/sky = 0.48) mask defined 
in Planck Collaboration XV (2013). The normalization is such 
that J df = /sky. The median integration time is about 0.42 s 
whereas the mean is 0.56 s. Notice that about 0.5 % of the pixels 
are observed for more than 5 seconds; these mostly lie around 
the Ecliptic poles. 



6.5. CO correction 

Some emission from the rotational transition lines of CO is 
present in all HFI bands. It is especially significant in the 100, 
217, and 353 GHz channels, due to the 115, 230, and 345 GHz 
CO transitions. This emission comes largely from the Galactic 
interstellar medium, mainly located at low and intermediate 
Galactic latitudes (\b\ < 20°). Due to the wide spectral cov- 
erage of Planck (Planck Collaboration IX 2013), we are able 
to produce velocity-integrated CO line maps for the first three 
transistions /=1^0, 2^1, 3^2. A complete description of the 
methods used, the validation tests, and the CO products are 
given in Planck Collaboration XIII (2013). These all-sky CO 
map products can be used for astrophysical studies of the in- 
terstellar medium or in component separation methods. 

Given the fraction of the sky used for the estimation of cos- 
mological parameters (/ky ^0.5) and that CO is only signific- 
ant in less than 1 % of the sky for these high Galactic latitudes, 
simply masking the CO-emitting regions is the course of action 
taken in cosmological studies with Planck. 

6.6. Spectral bands and colour corrections 

The HFI detector spectral response was determined through a 
series of ground-based tests conducted with the actual HFI focal 
plane in a simulated space environment, i.e., evacuated and cryo- 
genic, prior to launch. During these tests, the HFI focal plane 
was coupled to a broadband radiation source observed through 
a continuously scanned polarizing Fourier transform spectro- 
meter (FTS). The individual HFI detector responses to the mod- 
ulated FTS signal, i.e., the observed interferograms, were recor- 
ded alongside that of a reference bolometer located inside an 
integrating sphere within the experimental setup. Multiple scans 
were averaged together to increase the signal to noise ratio of 
the resultant spectra and to allow an estimate of the spectral 
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Figure 22. Power spectra of the diff'erence between intensity 
maps reconstructed from the diff'erence between the first and 
second half of each ring (denoted "HR dif"), or from the dif- 
ference between survey one and survey two. 



uncertainty. Standard FTS data processing techniques were em- 
ployed to obtain the observed spectra of each HFI detector and 
that of the reference bolometer; the ratios of the detector spec- 
tra against the reference bolometer spectrum provide the relative 
response of each HFI detector. This method incorporates the de- 
tector throughput (i.e., etendue AQ, see Planck Collaboration IX 
2013) in the spectral response as the reference bolometer ac- 
cepts a much broader field of view, i.e., approximately 27rsr 
over all frequencies. The HFI detector spectral transmission 
data, and their frequency band-average transmission spectra, 
are described more fully in Planck Collaboration IX (2013) 
and Planck Collaboration ES (2013), where a description of the 
band-averaging scheme is also provided. 

As discussed in Sect. 6.2, HFI is calibrated by the CMB 
dipole for the four 100-353 GHz bands, expressed in units of 
KcMB, and by planets for the two 545-857 GHz bands, ex- 
pressed in units of Wm"^sr"^Hz"^ (equivalently, MJysr"^). 
This unit conversion is defined to translate data between units 
of KcMB and MJy sr"^ , the latter for spectra with a profile con- 
sistent with a spectral index of -1, following the IRAS conven- 
tion (Beichman et al. 1988). Colour-correction coefficients are 
provided for a variety of spectral indices (see Fig. 25) and other 
common spectral profiles including modified blackbody profiles 
and planetary spectra. It is thus important to ensure that data pro- 
cessing, e.g., component separation, is conducted in a consistent 
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Figure 23. The zodiacal light contribution in HFI channels, with 100 GHz to 857 GHz from top to bottom. Units are yuKcMB for 
100-353 GHz and MJy sr"^ for 545 and 857 GHz. Left column: Estimate of the correction. Right column: The survey one-survey 
two difference map after correction, which can be compared with the bottom row of Figs. 15 to 20. While there is no visually 
perceptible difference in the CMB channels, note the visible cleaning of the 353-857 GHz maps. 
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Figure 24. Angular power spectrum of the zodiacal light contri- 
bution from each channel, as labeled. Even multipoles are shown 
with filled boxes, odd multiples with crosses. The systematic dif- 
ference reflects the approximate North/South symmetry of the 
zodiacal contribution. Spectra of other channels are given in 
Planck CoUaboration XIV (2013). 
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Figure 25. Band-averaged colour correction for various power- 
law frequency spectra, normalized to spectral index -1. The 
dashed lines provide estimates of the uncertainty of these cor- 
rections. 

manner in units appropriate for the frequency band. The HFI 
policy is to operate in units of MJy sr"^ , as this is valid for all of 
the HFI frequency bands whereas Kcmb is only valid for 100- 
353 GHz. The Kcmb to MJy sr"^ unit conversion coefficients are 
therefore given by the following ratios 
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where riy) is the relevant spectral transmission curve, b'^ = 
dBy/dT\TcMB is the Planck function derivative, and is the nom- 
inal band centre frequency. 

Other unit conversions are also provided, to ysz (the 
Sunyaev-Zeldovich Compton parameter, discussed in 
Planck Collaboration XXI 2013), Kkms"^ (CO transmis- 



sion; see Planck Collaboration XIII 2013), and Krj (brightness 
temperature). It is important to note that brightness temperature 
is not restricted to a Rayleigh- Jeans spectral profile, but is 
a convention for an alternate expression of an arbitrary flux 
density in temperature units (Rybicki & Lightman 1986). 

The HFI unit conversion and colour correction philosophy 
expresses equivalent flux density at a predefined nominal fre- 
quency depending on the spectral profile in question, rather than 
maintaining a set signal amplitude at a varying effective fre- 
quency, which would also depend on the spectral profile in ques- 
tion. This is primarily done to ease the requirements on compon- 
ent separation, as having all components expressed at the same 
effective frequency for all detectors within a frequency band is 
much preferred over dealing with a variety of different effective 
frequencies for each detector, and each component. 

More details on the unit conversion factors are provided 
in Planck Collaboration IX (2013), while tables of conver- 
sion factors for Planck detectors (only HFI channels for 
some spectral shapes, e.g., CO lines) are provided in the 
Planck Collaboration ES (2013). Figure 26 illustrates the band- 
averaged spectra for each of the HFI bands; similar, and re- 
lated, plots for individual HFI detectors are provided in the 
Planck CoUaboration ES (2013). 

A series of IDL scripts, named the UcCC routines, is avail- 
able to allow users of Planck data to determine unit conver- 
sion and colour correction coefficients, and uncertainties, for 
user-specified profiles in addition to those nominally provided 
in the public IMG. Further details on these scripts are provided in 
Planck Collaboration ES (2013). 

6.7. Polarization 

Together with the intensity maps, we derive polarization 
maps for each HFI frequency from 100 GHz to 353 GHz (cf. 
Planck Collaboration VIII 2013). The Q and U Stokes parameter 
maps show impressive signal dominated by Galactic dust emis- 
sion and CMB with noise of very low level and nearly white 
spectrum. Nevertheless, at the present state of the reconstruction, 
polarization maps are dominated at large scales by systematic ef- 
fects. As a consequence, they are not included in the present re- 
lease. However, these data will be used in forthcoming analyses 
of foreground polarization, and hence we describe the major sys- 
tematic effects remaining in the polarization data. 

The Planck scanning strategy does not allow recovery of the 
polarization signal for each detector independently. We use a 
combination of at least four detectors oriented at 45° in order 
to reconstruct the Q and U Stokes parameters. Any miscalibra- 
tion between the datasets will thus lead to a mix between Stokes 
parameters and, in particular, a leakage from intensity to polar- 
ization. As discussed in Planck Collaboration VIII (2013), abso- 
lute calibration is derived from the Solar dipole and gain vari- 
ations due to ADC non-linearities are corrected to a few tenths 
of a percent in the time domain. Final inter-calibration uncertain- 
ties between HFI bolometers are of the order of 0.4 % or less for 
the three lowest HFI frequencies (0.39, 0.28, and 0.21 % respect- 
ively at 100, 143, and 217 GHz) and rise to 1.35 % at 353 GHz. 
The induced leakage is dominated by the Solar dipole intens- 
ity, which produces large scale features in the polarization maps. 
We estimated the leakage due to calibration mismatch using sim- 
ulations. Figure 27 shows that the leakage is dominant over the 
CMB signal up to ^ = 100 in the EE angular power spectrum as 
predicted by the best-fit cosmological parameters measured from 
Planck temperature data (Planck Collaboration XVI 2013). 
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Figure 26. The band-averaged transmission spectra for all HFI frequency bands. The locations of the various relevant CO transitions 
are marked with vertical lines. 



In addition to the calibration mismatch, the difference 
in the bolometer spectral transmissions will also introduce 
a leakage from the dust emission intensity into polarization. 
Using the spectral transmission as measured on the ground 
(Planck Collaboration IX 2013), we estimated the level of leak- 
age introduced by bandpass mismatch as shown in Fig. 27. 



7. Data consistency and validation 



7.1. Consistency with a severe selection of data 



Polarization angles and efficiency have been assessed during 
ground measurements. Values together with their uncertainties 
are reported in Rosset et al. (2010). We derive uncertainties on 
angular power spectra using simulations for both polarization 
angle and efficiency independently. 



Figure 27 shows the level of specific sources of error in EE 
polarization angular power spectra compared to the CMB and 
the noise level (estimated from half-ring diff'erences) computed 
on 90 % of the sky excluding the Galactic plane. In the polar- 
ized HFI bands (100, 143, 217, and 353 GHz) systematic er- 
rors dominate over the CMB signal at low multipoles (£ < 100) 
due to miscalibration leakage and bandpass mismatch leakage. 
Polarization angle and polarization efficiency uncertainties are 
second order eff'ects. For Galactic emission polarization studies, 
the signal is more than one order of magnitude more intense. At 
higher multipoles, HFI polarization maps are noise-dominated 
(at £ > 200) but there is no evidence for correlated noise from 
channel to channel. 



Another check based on real data consists of comparison with 
data in which the pointing periods have been much more 
severely discarded. As mentioned in Sect. 3.10, the discarded 
pointing periods represent less than 1 % of the total integration 
time for the full mission and 0.2 % for the nominal mission, the 
diff'erence being due to the numerous solar flare events during 
the fifth survey. For this consistency test, we discard more than 
30 % of the data of the full mission, following much more severe 
criteria. For example, we use very conservative thresholds on the 
stationarity of the noise (details in Appendix B). We then create a 
map for each frequency with this selection. We compute pseudo- 
spectra of both the normal and the severely discarded selection. 
As can be seen in Fig. 28, the diff'erence is compatible with the 
Gaussian case for £ > 500. 



There is a rise in the pseudo- spectrum for ^ < 500. However, 
such a feature is far below the cosmic variance (at the percent 
level). So the level of non-stationarity of the noise in the data 
of the nominal mission is negligible and does not impact the 
cosmological results. 
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Figure 27. Uncertainties on polarized power spectra due to residual systematics in the HFI polarization maps compared to the EE 
spectrum predicted for the best-fit model from Planck temperature data. 
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Figure 28. The diff'erence of the pseudo- spectra for the normal 
and the severe selections, in units of the expected standard de- 
viation. If our assumptions are correct, the residual noise is sup- 
posed to have a Gaussian distribution centred on zero with a dis- 
persion of one (dashed lines). The red points correspond to a 
statistical check with randomly discarded rings. The points are 
band-averaged in bins of width ~ 91. 



7.2. Difference map consistency tests 

We perform a series of diff'erence map consistency tests to eval- 
uate the contribution of residual systematic eff'ects to the angu- 
lar power spectrum of the temperature anisotropics. We gener- 
ate diff'erence maps of halves of the nominal mission data in 
which the signal is expected to fully subtract out in the ab- 



sence of systematic eff'ects. The resulting maps are then subtrac- 
ted from systematics-free Yardstick realizations (Appendix A) 
of the same diff'erences propagated through the same mapmak- 
ing pipeline as that used for data processing. The amplitude of 
the power of the remaining signal in each likelihood bin can 
thus be used as an estimate of the level of the residual sys- 
tematic contamination in this bin. This estimate is derived from 
cross spectra between diff'erent detector sets to allow for a direct 
comparison to the spectra used as inputs by the likelihood code 
(Planck Collaboration XV 2013). The Yardstick simulations are 
described in detail in Appendix A. 

Many diff'erences of halves of the data can be expected to be 
compatible with an absence of signal. Given Planck's scanning 
strategy, we expect the bulk of the remaining systematic con- 
tamination to be captured in survey diff'erence maps. In partic- 
ular, transfer function errors, far-sidelobe pickup, gain instabil- 
ity, pointing drifts, and residual glitches will leave an imprint 
in these maps. Half-focal plane diff'erence maps will be par- 
ticularly sensitive to relative gain errors, while also providing 
some sensitivity to sidelobe pickup. Finally, diff'erence maps ob- 
tained by subtracting PSB-only from SWB-only sky maps probe 
any unexpected behavioural diff'erence (e.g., beam mismatch) 
between polarization- sensitive bolometers and unpolarized dec- 
tectors, and provide a check of the (/, Q, U) decomposition ap- 
plied to the data from the former. Bandpass diff'erence mises- 
timation between the two detector technologies will also leave a 
signal in these maps. 

Since the survey diff'erence test is the most stringent diff'er- 
ence map consistency test at our disposal, we focus the discus- 
sion in this section on the corresponding results. We did not de- 
tect null test failures in any of our other tests from the list above. 
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Figure 29. (Survey 2 - Survey l)/2 consistency test result from the dsl x ds2 cross spectra at 100, 143, and 217 GHz with the 
/sky - 0.30 mask derived from that used in the primary cosmological analysis (see the text). Left: the residual signal in each 
bandpower up to the highest multipole used in the likelihood code (Planck Collaboration XV 2013). For clarity, the bins used in the 
likelihood analysis are grouped four-by-four for £ > 60 and eight-by-eight for £ > 1250 in the plots above. We show in each bin 
the variance of the Q distribution from the simulations. Note that the amplitude of the residual signal is under the binned sample 
variance envelope expected for the /sky - 0.30 mask (shown in black) up to ^ ~ 1000 at all three frequencies. Right: zoom-in on the 
low-^ part of the spectrum, showing that although the difference test failures are highly significant at these scales, their amplitudes 
at all three CMB frequencies are less than or of the order of 0.5 fiK^ for all multipoles up to ^ = 200, which is a tiny fraction of the 
binned sample variance at these scales. 



The nominal mission survey diff'erence test results are shown in 
Fig. 29, while Fig. 30 shows the outcome of a suite of diff'erence 
map consistency tests at 217 GHz. 

In Fig. 29, the angular power spectrum of the residual signal 
is obtained by taking the cross spectrum between detector sets 
dsl and ds2 at each frequency (see Table 3) after application 
of the /sky = 0.48 mask used in the primary cosmological ana- 
lysis (Planck Collaboration XV 2013) combined with a mask of 
the regions not seen by either the first or the second survey and 
a mask of the point sources included in the input PSM maps. 
The resulting sky fraction is 0.30 at 100 GHz, and 0.29 at both 
143 GHz and 217 GHz. Results are then binned according to the 
likelihood binning. For clarity, bins at multipoles higher than 
£ = 60 are grouped four-by-four, whereas those at ^ > 1250 
are grouped eight-by-eight. The plot at left shows the residual 
signal over the whole range of multipoles used in the likelihood 
analysis. At right, the same results are shown with a low-^ zoom 
(£ < 200). In both cases, we report the variance of the Q distri- 
bution in each bin as computed from the simulations. 

At all three CMB frequencies, the £ < 200 residuals are strik- 
ingly small, never exceeding about 0.5 yuK^. Excluding a single 
bin in each case, residuals actually remain below or at the level 
of 0.2 yuK^ over this range of multipoles. Although these non- 
zero diff'erences are detected with very high statistical signific- 
ance, they are many orders of magnitude smaller than the binned 
sample variance (shown in black in Fig. 29) at these scales. 
In particular, they cannot aff'ect the cosmological analysis. This 
stays true all the way up to ^ ~ 1000, at which point residuals 
become higher than binned sample variance at 100 and 217 GHz. 
At 143 GHz, this does not happen until £ reaches 1500. 

In the multipole range from 1000 to 2500, although the amp- 
litude of the residuals is almost always greater than the binned 
sample variance, the variance in the simulation results is signi- 
ficantly larger than in the £ < 1000 regime. As a result, the 100 
and 143 GHz residuals for that range of multipoles are fully com- 



patible with zero. However, the 217 GHz residuals in the same 
multipole range are not, as can be directly inferred from the plot, 
where a significant oscillatory feature starts at ^ = 1000. 

This survey-diff'erence test has been performed at each fre- 
quency for all combinations of the input maps used in the like- 
lihood analysis (two at 100 GHz, five at 143 GHz, and six at 
217 GHz), and for two survey diff'erences: Survey 1 - Survey 2 
and Survey 1 - Survey 3. In addition to the 217-dsl x 217-ds2 
cross spectrum shown in Fig. 29, only two other cross-spectra 
fail this test, namely 217-1 x 217-ds2 and 217-1 x 217-3. In the 
likelihood analysis (Planck Collaboration XV 2013), we have 
checked that the determination of the cosmological parameters is 
not aff'ected by the inclusion, or not, of these three cross-spectra. 

7.2.1. ADC non-linearity impact 

We check the impact of including the eff'ect of ADC non- 
linearity by the comparison of data and simulations of the two 
PSB pairs at 143 GHz. Note that this eff'ect had been suspected 
to be significant, but was neither corrected for in the nominal 
mission products, nor included in the corresponding Yardstick 
simulation, since the relevant modelling information was only 
obtained after the end of the HFI cryogenic phase through dedic- 
ated data gathering during the warm phase. Figure 3 1 shows that 
the inclusion of this eff'ect leads to similar ^ < 100 deviations 
from the Yardstick compared to those observed in the survey 
diff'erence consistency test of the 143 GHz data. The inclusion 
of this eff'ect makes little diff'erence at higher values of £. 



7.2.2. Noise estimation bias from half-ring maps 

A detailed comparison of the power spectra of the sum and dif- 
ference of half-ring maps shows a constant off'set of a few per- 
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Figure 30. Examples of map difference consistency tests at 
217 GHz with the /sky - 0.30 mask derived from that used in 
the primary cosmological analysis (see the text). The range of 
multipoles shown and the binning used are identical to those de- 
scribed in Fig. 29. We again show in each bin the variance of 
the Q distribution from the simulations, and in black the binned 
sample variance for our mask. "Left" and "Right" refer to the 
two halves of the focal plane, whereas "In" (respectively, "Out") 
refer to detectors toward the centre (respectively, the periphery) 
of the focal plane. Unlike in the survey difference test shown in 
Fig. 29, no significant null test failure is detected in these tests. 
The detector combinations used in the example plots above are 
(5 - 8) X (1 - 4), (7 - 8) X (3 - 4), and (dsl - 3) x (ds2 - 4) for, re- 
spectively, the (left-right), (in-out), and (PSB-SWB) tests. As 
in Fig. 29, these power spectra are divided by four to allow for a 
direct comparison of the amplitude of the residuals to that of the 
components in the nominal mission power spectrum. 



cent at large multipole. At these scales, there is no longer any 
signal remaining in the sum map, due to the beam cut-off (see 
Fig. 32 for an example), and we would expect that the difference 
map spectrum would give a precise determination of the high-^ 
tail of the sum map, while it is actually slightly below. This off- 
set shows that a noise estimate from half-ring maps difference is 
a slightly biased indicator of noise in the sum maps. This results 
from the deglitching algorithm. 

In the first step, the pipeline estimates the signal by averaging 
the timelines into rings. This signal estimate is removed in order 
to aid detection of glitches by comparison with the noise level. 
Each positive sample at least 3cr above the noise level is flagged 
as a glitch. One side effect of this thresholding is that the pos- 
itive tail of the timeline noise distribution is now clipped at 3cr. 
Note that the signal estimate within a ring is of course still noisy, 
so that its removal from the timeline before deglitching subtracts 
this residual noise component periodically (coherently at all har- 
monics of the spin frequency), in effect periodically modulating 
(in both half-rings) the flagging and clipping level of the noise 
distribution. In the difference of half-ring maps, this common 
contamination disappears, while it is present in the sum map, 
leading to this slight bias. 

Figure 32 shows that the inclusion of the deglitching step 
in simulations accounts quantitatively for the difference in the 
power spectra of the sum and difference of half-ring maps. By fit- 
ting the ratio of the spectra of sum and difference maps (masked 
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Figure 31. Simulated effect of the ADC non-linearity on S1-S2 
survey difference at 143 GHz. Top: Difference test of a single 
realization of the full pipeline, assuming a linear response of the 
ADC. Middle: the same, but with our current best model of the 
non-linearity response of the ADC. This can be compared with 
the bottom panel, which shows the result of the test applied to the 
delivered data in exactly the same conditions. This demonstrates 
that indeed the effect of ADC non-linearity can explain most of 
coherent deviation seen at ^ < 100. 



to retain 37 % of the sky) at very high £ (around £ ~ 7000) one 
finds the level of the bias which are given in table 2. 
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Figure 32. Comparison of the high-^ part of the power spectra 
of the sum and difference of half-ring maps for the 143-la de- 
tector, both for flight data and simulations where glitches were 
introduced and deglitched with the actual processing pipeline. 
The comparison shows that the simulations provide a good de- 
scription of this effect (which is not included in the Yardstick 
simulations). Top\ the sum and difference spectra. Bottom', ratios 
of the sum and difference spectra for both data and simulations, 
which differ by about 2 %. 



Table 2. Percentage increase of the power spectrum of the dif- 
ference of half-ring maps needed to improve the accuracy of the 
noise level estimate in frequency maps. 
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0.61 + 0.08 % 0.75 + 0.09 % 
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7.3. Power spectrum consistency tests 

A powerful check of the data follows from assessing the 
signal consistency between detectors at the same frequency 
given the prior information on the effective beam uncertainties 
and calibration. At each frequency, an estimate of an adapted 
power spectra model is processed with the SMICA likelihood 
(Planck Collaboration XV 2013). At 143 and 217 GHz where the 
CMB is dominant, the spectrum on the measured sky is mod- 
elled as aa^Q + n^, a sum of the the free signal power, Q, plus 
a free independent noise spectrum, n^, where the known vec- 
tor a accounts for the fixed CMB emissivity at each frequency 
(bold-faced vectors have elements at each frequency). At higher 
frequencies, dust emission begins to dominate the CMB, and 
some foreground modelling is necessary. Thus the CMB power 
is fixed to the angular spectrum determined at 143 GHz, the noise 
is kept free and a free foreground is added, with colour cor- 
rection determined from a grey body spectral shape integrated 
though the measured spectral response of each detector [40 % of 
the sky is used to estimate the binned empirical power spectrum 
with spice (Szapudi et al. 2001)]. Furthermore, beam errors are 
modelled as a linear combination of templates, thanks to a beam 
error mode analysis (Planck Collaboration VII 2013). Then, five 



(4) 



where is the modelled spectrum, the free beam model, 
b the free dust emmisivity at each frequency, g£ the free dust 
power, and the free noise power. The beam model for detector 
d is written Bf[d] = exp(2^ ^'^7^^'^) where 6 are the beam- 
eigenmode amplitude parameters and T the templates. Due to 
the intrinsic ambiguity between spectral power and beams, beam 
simulations are used to determine the variances of a Gaussian 
prior on these parameters. Finally, at 100 GHz, with only two 
detector sets, three parameters per £ (one signal and two noise 
power) allows an estimate of the empirical spectrum, but relative 
calibration and beam parameters estimation are then unfeasible. 

Figure 33 shows the result of such an analysis. At 143 and 
217 GHz, where the CMB is dominant and used as calibrator, 
the only adjustments necessary are a gain recalibration by a very 
small amount (0.2 % at most) and the addition of a fraction of a 
(T of the first beam uncertainty eigenmode. At 545 and 857 GHz, 
where the CMB is subdominant and planets are used as calib- 
rators, the required recalibration remains just as small (about 
0.2 %), but larger beam adjustments are necessary. The 353 GHz 
calibration is the one requiring the largest changes, of about 1 %. 
At 217 GHz and to a lesser extent 353 GHz, the SWB and PSB 
beam window functions are at opposite sides of the "optimal" 
one. In conclusion, by properly modelling foregrounds at high 
frequency, cross-checks on signal anisotropy with 40 % of the 
sky shows that map are well calibrated (relatively) within a fre- 
quency. 

Figure 34 shows the magnitude of the beam corrections de- 
termined by maximising the internal consistency for the most 
sensitive CMB channels of Planck. The largest corrections 
would be about around ^ ~ 100 and about around 

the third peak of the CMB spectrum (£ ~ 800). For the delivered 
frequency maps, which combine all detectors at the same fre- 
quency, the effect is even weaker, and for the likelihood analysis, 
this is fully accounted for (these beam eigenvalues are either ex- 
plored or marginalised over, including in addition their expected 
cross-correlations). 

7.4. Consistency of the calibration across frequencies 

In the previous section, we looked at the consistency of de- 
tector spectra within a given frequency band. A similar ana- 
lysis of the spectra used in the high-^ part of Planck likelihood 
Planck Collaboration XV (2013) shows that the best relative re- 
calibration (for the mask considered and accounting for the first 
beam eigenmode), between all 13 detectors maps in the 100, 143 
and 217 GHz channels, would be less than 0.2%. This level has 
no sizeable effect on the constraints on cosmological parameters. 

Finally, we can further broaden our consistency checks 
by using the SMICA component separation method from 
44 to 353 GHz, in effect intercalibrating on the common 
CMB anisotropics themselves, taking 143 GHz as the ref- 
erence (see Sect. 7.3.3 of Planck Collaboration VIII 2013 
and Planck Collaboration XII 2013). Figure 35 compares 
the results of such an analysis^ on an FFP6 simulation 
(Planck Collaboration ES 2013) and on the 2013 maps, using a 



^ The analysis is actually done over 40 % of the sky, for various £ 
ranges; the results shown here show the average in the multipole range 
[100, 679] for the simulations. 
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Figure 33. Calibration and beam uncertainty mode determination for all detectors (or detector sets) within a frequency. The grey 
areas show thelcr priors on each beam parameters (without assuming a prior for gain calibration, The first detector in each series is 
used as a fixed reference). The coloured arrows show thelcr posteriors. The bottom right shows the variations of the power spectra 
in the primary CMB channels when applying the recalibration and beam uncertainty correction described in the text. Altogether, this 
leads to variations smaller than 0.4 % at all £ < 2500. Note that the shape variations with £ are very much smaller (A < 5 x 10""^). 



five-component foreground model. The simulations accurately 
validate the approach, in a realistic setting (albeit with decreas- 
ing accuracy at 44 and 353 GHz). For flight data, the multipole 
range over which the calibration factor is averaged has been 



split in two, a relatively low-^ range (50-300), more appropri- 
ate to include the noisier 44 and 70 GHz channels, and a higher- 
£ range, (300-700), more appropriate for HFI data. The figure 
confirms that the 100-217 GHz map calibrations are consistent 
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Figure 34. This shows the derived overall beam correction of 
each detector set which minimizes their mismatch, given the 
model used, at 143 and 217 GHz where there are sufficient data 
to simultaneously estimate beam parameters along with sig- 
nal and noise power spectra. These corrections peak around 
i ~ 200 and 800 and are at most of the order of three yuK^. 



to 0.2 %, as discussed earlier. The 353 GHz results would sug- 
gest a slightly higher calibration than the one we finally adopted 
(but the SMICA analysis does not account for the Rayleigh cor- 
rection (see Planck Collaboration VIII 2013) and we may also 
reach the limit of a five-component foreground model). This ana- 
lysis also suggests a residual calibration diff'erence between LFI 
at 70 GHz and HFI at 100 GHz of roughly 0.3 %. 

Further checks involving comparison with the LFI and 
WMAP data are discussed in Planck Collaboration XI (2013). 

7.5. Further validation 

Other tests are described in various other papers: 

- Planck Collaboration XI (2013) compares the maps from the 
HFI and LFI instruments on Planck and from the WMAP 
satellite. We demonstrate excellent consistency amongst the 
three sets of maps produced with diff'erent instruments over a 
a factor of thirty in wavelength, and especially over the nar- 
row range of 70-100 GHz covered by all three instruments. 
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Figure 35. Recalibration factor maximising the CMB consist- 
ency in in simulations (black) and in the data considering diff'er- 
ent multipole ranges (red and blue), at each Planck frequency in 
GHz given on the horizontal axis. 



illustrating the excellent control of systematic errors we have 
achieved with Planck. 

- Planck Collaboration XII (2013) discusses diff'erent meth- 
ods for separating the various astrophysical and cosmolo- 
gical data. Methods with very diff'erent assumptions about 
the physics of the diff'erent components, and about the nature 
of the noise in the Planck maps, are able to converge on equi- 
valent maps for diff'erent components, most importantly of 
the CMB itself. 

- Planck Collaboration XV (2013) discusses the angular 
power spectrum, Q, computed from Planck data, as well as 
the likelihood function that embodies its full probability dis- 
tribution. The likelihood is computed by taking into account 
explicit models foreground contamination and of systematic 
eff'ects such as beam errors and the relative calibration of 
diff'erent detectors. The resultant spectra explicitly illustrate 
consistency across those bands with the highest signal-to- 
noise for the cosmological spectra at 100, 143, and 217 GHz, 
at the level of power spectra and of cosmological parameters. 
As discussed as well in 7.2, the likelihood is robust to the in- 
clusion or otherwise of even the 217 GHz detectors that show 
some mild non-zero features in the spectra of their diff'erence 
maps. 



8. Summary of product characteristics and 
conclusions 

The data provided to the Planck collaboration by the HFI DPC, 
in addition to the I MO, consisted of the following (letters refers 
to entries in the final column of Table 3): 

1. C: beam-corrected power spectra for each CMB detector set 
used in the Planck likelihood code, together with a descrip- 
tion of the relevant beam uncertainties; 

2. N: frequency maps for the nominal mission duration, not cor- 
rected for ZLE/FSL; 

3. B: eff'ective beams and uncertainties (for N); 

4. F: spectral bands and uncertainties; 

5. H: hit count maps (for N); 
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Table 3. Detector sets used in this data release. Unpolarized SWBs are used alone, while PSBs (with individual bolometers denoted 
"a" and "b") are used by pair. Here we consider only the temperature map extracted from the analysis of four bolometers (two 
pairs of PSBs). The weights indicate the relative weighting used in producing maps out of the TOI of several detectors. The relative 
weights in a set are given in the numeric order of each detector (e.g., la, lb, 2a, 2b, 3a, 3b. . . for the 100-dsO set). The last column 
details the specific products created for each set. "N" stands for maps created from nominal mission data (with their beam and band 
information, "B" and "F"), "H" is the corresponding hit map per pixel, "SI" and "S2" stands for the survey maps generated from 
the data collected during the first six month or during the next six months. "HRl" and "HR2" are for maps made out from the first 
or second half of each ring. "C" stands for the detector-set spectra, corrected for the beam transfer function, for the detector sets 
used in the high-^ likelihood. "Z" corresponds to the ZLE/FSL corrected frequency maps. 



Set name 


Frequency 


Type 


Detectors in the set 


Weights in the set 


Products 




[GHz] 










100-dsO 


100 


MIX 


All 8 detectors 


0.33,0.44, 1.50, 0.76 


N, B, F, H, SI, S2, HRl, HR2 










2.08, 1.40, 0.97, 30.52 


Z-N, Z-S1,Z-S2 


100-dsl 


100 


PSB 


la+lb + 4a+4b 


0.58, 0.78, 1.71,0.93 


C 


100-ds2 


100 


PSB 


2a+2b + 3a+3b 


1.05,0.53, 1.45,0.98 


C 


143-dsO 


143 


MIX 


11 detectors 


1.04, 0.44, 1.06, 0.93,0.88,0.94 


N,B,F,H, S1,S2, HR1,HR2 










0.80, 0.66, 1.32, 1.27, 1.66 


Z-N, Z-S1,Z-S2 


143-dsl 


143 


PSB 


la+lb + 3a+3b 


1.26, 0.53, 1.06, 1.15 


c 


143-ds2 


143 


PSB 


2a+2b + 4a+4b 


1.23, 1.08, 0.92, 0.76 


c 


143-ds3 


143 


SWB 


143-5 




c 


143-ds4 


143 


SWB 


143-6 




C 


143-ds5 


143 


SWB 


143-7 




C 


217-dsO 


217 


MIX 


All 12 detectors 


1.54, 1.44, 1.62, 1.83,0.60, 0.77 


N,B,F,H, S1,S2, HR1,HR2 










0.62, 0.67, 0.81,0.85,0.64, 0.59 


Z-N, Z-S1,Z-S2 


217-dsl 


217 


PSB 


5a+5b + 7a+7b 


0.79, 1.02, 1.07, 1.12 


C 


217-ds2 


217 


PSB 


6a+6b + 8a+8b 


0.98, 1.06, 1.01,0.94 


C 


217-ds3 


217 


SWB 


217-1 




C 


217-ds4 


217 


SWB 


217-2 




C 


217-ds5 


217 


SWB 


217-3 




C 


217-ds6 


217 


SWB 


217-4 




C 


353-dsO 


353 


MIX 


All 12 detectors 


2.45,2.38, 0.44, 0.65,0.61,0.57 


N,B,F,H, S1,S2, HR1,HR2 










0.64, 0.64, 0.31,0.34, 1.62, 1.35 


Z-N, Z-S1,Z-S2 


353-dsl 


353 


PSB 


3a+3b + 5a+5b 


0.75, 1.09, 1.09, 1.07 




353-ds2 


353 


PSB 


4a+4b + 6a+6b 


1.33, 1.25,0.68,0.73 




353-ds3 


353 


SWB 


353-1 






353-ds4 


353 


SWB 


353-2 






353-ds5 


353 


SWB 


353-7 






353-ds6 


353 


SWB 


353-8 






545-dsO 


545 


SWB 


3 detectors (1,2, 4) 


0.94, 1.10, 0.96 


N,B,F,H,S1,S2, HR1,HR2 












Z-N, Z-S1,Z-S2 


857-dsO 


857 


SWB 


All 4 detectors 


1.14, 1.11, 1.15,0.60 


N,B,F,H, S1,S2, HR1,HR2 



Z-N, Z-S1,Z-S2 



6. HR: half-ring maps, HRl and HR2, two maps made from 
the first and the second half of each pointing period of the 
nominal mission; 

7. S: survey maps, SI and S2, made from the data of survey one 
only and survey two only; and 

8. Z: zodiacal light emission and far sidelobe corrected fre- 
quency maps and their survey maps (Z-N, Z-Sl, Z-S2). 

Table 3 summarize the various detectors sets from which we 
made these products while Table 4 summarizes the map char- 
acteristics needed to make use of the maps. Further details and 
comments are given below. 

Angular Response 

Eff'ective Beams provide the response of a map pixel to the sky 
(lines hl-b7 of Table 4). They are based on a determination of 
the scanning beam from planet scans (Sect. 5.1), which include 



the eff'ect of the optical beam, the electronic detection chain, and 
the TOI processing pipeline. Scanning beams are an important 
intermediate product, although they are usually not relevant for 
astrophysical applications. The eff'ective beam further accounts 
for the combined eff'ect of the scanning strategy and additional 
data processing. 

Different applications need diff'erent levels of accuracy and 
detail, from the mean full-width at half-maximum (FWHM) of 
a symmetrical Gaussian description (line b4) to an actual point 
spread function at each specific map location. Beam solid angles 
are given in line bl . Mean ellipticities are given in line b6. Lines 
b2, bS, and b7 provide information on the angular response un- 
certainties and variation across the sky. Our complete uncer- 
tainty budget, including covariances of the beam eigenvalues 
between detectors, in a form which is usable for the power spec- 
tra analyses of Planck Collaboration XV (2013), is provided in 
Planck Collaboration VII (2013). 
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Table 4. HFI Nominal Maps — Main Characteristics. 



Quantity 


Units 














Notes 


Reference frequency y 


[GHz] 


100 


143 


217 


353 


545 


857 


a\ 


Number of bolometers 








1 1 
1 1 


iz 


iz 


a 
J 


A 


al 


Effective beam solid angle Q 


[arcmin^] 


105.78 


59.95 


28.45 


26.71 


26.53 


24.24 


bl 


Error in solid angle ctq 


[arcmin^] 


0.55 


0.08 


0.03 


0.02 


0.03 


0.03 


bl 


J.' t ' J.' / \ A 

Spatial variation {rms) Ail 


r • 2i 

[arcmm^J 


0.31 


0.25 


0.27 


0.25 


0.34 


0.19 


b3 


Effective beam FWHMi 


[arcmin] 


9.66 


7.27 


5.01 


4.86 


4.84 


4.63 


b4 


Effective beam FWHM2 


[arcmin] 


9.65 


7.25 


4.99 


4.82 


4.68 


4.33 


b5 


Effective beam ellipticity € 




1.186 


1.036 


1.177 


1.147 


1.161 


1.393 


b6 


Variation {rms) of the ellipticity A 6 








U.UJU 


U.Uzo 




U.U/D 


bl 


Sensitivity per beam solid angle 


[vt/K] 


10 


6 


12 


39 






c\ 




[kJy sr-i] 










13 


14 


cl 


Sensitivity 


[)t/Kdeg] 


1.8 


0.8 


1.0 


3.5 






c2 




[kJy sr~^ deg] 










1.1 


1.1 


cl 


Relative calibration accuracy 


[%] 


< 0.2 




< 0.2 


< 1 


< 5 


< 5 


d\ 


Absolute calibration accuracy 


[%] 


<0.5 


<0.5 


<0.5 


< 1-2 


< 10 


< 10 


dl 


Galactic zero-level offset 


[MJy sr-i] 


0.0047 


0.0136 


0.0384 


0.0885 


0.1065 


0.1470 


el 


Galactic zero-level uncertainty 


[MJy sr-i] 


0.0008 


0.0010 


0.0024 


0.0067 


0.0165 


0.0147 


el 


CIB monopole prediction 


[MJy sr-i] 


0.0030 


0.0079 


0.033 


0.13 


0.35 


0.64 


el 



Channel map reference frequency, and channel identifier. 
Number of bolometers whose data was used in producing the channel map. 
Mean value over detectors at the same frequency. See Sect. 5. 
As given by simulations. 

Variation {rms) of the solid angle across the sky. 

FWHM of the Gaussian whose solid angle is equivalent to that of the effective beams. 
FWHM of the mean best-fit Gaussian. 

Ratio of the major to minor axis of the best-fit Gaussian averaged over the full sky. 
Variability {rms) on the sky. 

Estimate of the noise per beam solid angle given in hi. See Sect. 6.3. 
Estimate of the noise scaled to one degree assuming the noise is white. 
Relative calibration accuracy between frequency channels. See Sect. 6. 
This estimate includes the estimated uncertainty of the calibrating source. 

Offset to remove at each frequency to set the Galactic zero level. The values quoted correpond to v/^ 
the conversion factor is about 244, 371, 483, 287 MJy sr'^ K'J^g. 
Overall error on the map zero point, for a y/y = const spectrum. 

According to the Bethermin et al. (2012) model, whose uncertainty is estimated to be at the 20 % level. (Also for vly 



const. At the four lowest frequencies. 



: const). 



Sensitivity 

The numbers given in line cl of the table indicate the rms contri- 
bution of the noise per beam solid angle, for the median integra- 
tion time, as scaled from maps at the full resolution (pixels of 1 .7 
arcmin). The numbers given in line c2 of the table also converts 
this number into an rms of the noise in pixels of one degree on a 
side, again assuming the rms varies inversely proportional to the 
pixel scale, as if it were white noise. While convenient to com- 
pare the sensitivity of different experiments, this number may be 
misleading, since it is actually an overestimate of the real noise 
in the map at the one-degree scale. 

It is interesting to note that the sensitivity of the 
map delivered is quite comparable to pre-flight predictions 
(Planck Collaboration 2005), despite our flagging out of a signi- 
ficant fraction of data. This decrease of integration time per pixel 
is approximately compensated by the better-than-required sens- 
itivity of the individual bolometers. A more elaborate descrip- 
tion of the noise is made available as "half-ring" maps whose 
differences offers a quite accurate view of the small scale noise, 
differint by up to half a percent. 



Photometric Accuracy 

The photometric calibration (Sect. 6.2) of the 100 to 353 GHz 
channels relies on the Solar dipole. Comparing the common 
CMB component at these frequencies shows that the relative ac- 
curacy between these channels is better than 0.2 % between 100 
and 217 GHZ, and at the percent level for the 353 GHz chan- 
nel. For the two highest frequencies, this relative accuracy is at 
the five percent level (lines dl-dl). We have also estimated the 
zero level offset of the maps which needs to be accounted for in 
Galactic and cosmic infrared background studies (lines el-e3). 

Spectral response and conversions 

The accuracy of the HFI spectral response characterization 
(Sect. 6.6) is validated using a variety of HFI in-flight obser- 
vations. Comparisons of zodiacal light observations (Sect. 6.4) 
with the bandpass data reveal an out-of-band signal rejection of 
better than 10^. The bandpass-based unit conversions and colour 
corrections were compared against those derived using sky-only 
data for sources including Sunyaev-Zeldovich (SZ) clusters, dust 
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emission, and CO emission. The SZ observations show agree- 
ment with the bandpass data within the quoted uncertainty, as 
do the 100, 217 and 353 GHz dust coefficient comparisons. The 
143 GHz dust coefficient comparison shows differences between 
the sky and bandpass calibration methods employed; however, 
the dust component is not dominant at 143 GHz so this dis- 
crepancy, although important to understand, is not alarming. 
The CO-based dust colour correction coefficients show good 
agreement with their bandpass-based counterparts, while there 
are some differences between the bandpass-based and sky-based 
CO unit conversion coefficients. These differences comprise 
areas of ongoing effort which will improve the understanding 
of the HFI instrument and data; this will in turn lead to better 
data analysis and improved knowledge of systematic and cal- 
ibration uncertainties. Details are provided elsewhere, primarily 
in Planck Collaboration IX (2013) and Planck Collaboration ES 
(2013). 

Data limitations, validation, and checks for systematic errors 

The most powerful top-down tests are the statistics of differences 
between sky maps. The power spectra of these difference maps 
have been computed in the same way as those used in the Planck 
CMB likelihood (which are based on the nominal maps). We 
have applied a large number of such tests, differencing maps 
made from various detector groups (in/out, left/right), and sky 
surveys. The sky survey differences maximize the ability to de- 
tect effects potentially affecting a specific detector or detector 
set; they have been applied to all likelihood inputs. We have 
found three potentially significant departures at large multipoles, 
for cross-spectra involving 217 GHz detectors. A jackknife test 
in the CMB power spectrum and likelihood paper verified that 
their removal does not affect significantly the cosmological con- 
straints (Planck Collaboration XV 2013). 

We also found small, microkelvin-level, residuals at low 
multipoles, which although very significant statistically, are in- 
significant for a cosmological analysis of the temperature. We 
have verified through simulations that most of this effect is due 
to the limitations of the partial correction made for the ADC 
non-linearity. This limits the accuracy of our dipole calibration 
to about 0.2 % in the CMB channels, and prevents us from using 
polarization information at low multipoles. However, the data 
collected at the end of the cryogenic chain will allow a more 
satisfactory correction for ensuing Planck data releases. 

Conclusions 

The six intensity maps and ancillary information presented 
above constitute an unprecedented source of information in this 
frequency range, which is in line with our pre-flight expecta- 
tions. The level of understanding reached while analysing this 
release bodes well for the expected next release of Planck data 
in 2014, which will take aim at polarization on all scales, and 
use nearly twice as much data. 
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Appendix A: Yardstick simulations 

The Yardstick simulation pipeline has been developed for two 
main purposes: to quantify the level of residual systematic ef- 
fects in the maps produced by the HFI DPC, and to validate our 
correction (or lack thereof) of the data for the systematic effects 
that we are able to model. The former is achieved by comparing 
combinations of the data from which the sky signal is expec- 
ted to vanish to an ensemble of signal and noise realizations of 
the same combinations of data subsets propagated through the 
same DPC processing as the data. A mismatch between the two, 
e.g., in angular power spectrum space, reveals a residual system- 
atic effect in the data. These tests are described in Sect. 7.2. In 
addition, validating the correction of a given systematic effect 
requires the implementation of both a model of the effect and its 
correction (the latter is of course unnecessary if the aim of a test 
is to justify an absence of correction). In Sect. 7.2.1, the valid- 
ation of the non correction in the data of the ADC nonlinearity 
is an example of such a test. The Yardstick pipeline has been 
designed in a modular way to allow for both kinds of tests. 
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For all simulations, the sky signal is taken from the Planck 
Sky Model (PSM), which has been described in detail in 
Delabrouille et al. (2012). This parametric model allows the 
generation of all-sky temperature and polarization maps of the 
CMB, the SZ effects, and diffuse Galactic emission (in partic- 
ular synchrotron, free-free, and thermal dust) with a resolution 
of a few arcminutes at all Planck frequencies. The PSM also 
includes an extensive point-source catalog, as well as spinning 
dust, CO line, and HII region models. From the point of view 
of data validation, a particularly useful feature of the PSM is its 
ability to generate random realizations constrained to match ob- 
servational data within their uncertainties. Although these realiz- 
ations cannot be fully independent (with the exception of those 
of the CMB), they do allow for the propagation of some sky 
modelling uncertainty through the DPC pipeline. The PSM also 
handles convolution with the spectral response of the HFI bolo- 
meters. 

The Yardstick pipeline uses the a^m spherical harmonic coef- 
ficients from the PSM as inputs. After generating the point- 
ing timelines directly from Planck's attitude history file, it 
feeds both to the LevelS-Core pipeline, software developed 
jointly by the HFI and LFI teams and documented in the HFI 
data processing early paper (Planck HFI Core Team 2011b). 
Levels -Core returns the beam-convolved signal TOI for each 
bolometer, which is fed back into the Yardstick pipeline. 
Noise TOI are generated by a dedicated module and include 
both white noise and 1// components. The high-frequency 
component is taken to be the detector noise described in 
Planck Collaboration VIII (2013), while the 1// noise is mod- 
elled following Planck HFI Core Team (2011a). The overall 
amplitude of the noise is given on a ring-by-ring basis by half- 
ring map null tests. The latter is further discussed in Sect. 7.2.2. 
The noise model does not include correlations between detect- 
ors. 

Instrumental effects can be added to the simulated data, each 
by the addition of a dedicated module to the pipeline. We have 
used this ability to improve our understanding of pointing and 
beam issues, as well as to test various aspects of glitch and 4 K 
line removal, among others. Section 7.2.1 shows an example of 
the use of the Yardstick simulation facility to evaluate the im- 
pact on null-test residuals of the non-correction of the ADC non- 
linearity. When evaluating the quality of the correction of a sys- 
tematic effect, we simply apply the relevant TOI processing prior 
to proceeding with mapmaking. The impact of the transfer func- 
tion convolution and deconvolution is in any case accounted for. 

Finally, TOI are propagated through the DPC mapmaking 
and calibration pipelines following the procedures described in 
this paper, including gain correction and application of the vari- 
ous flags (Planck Collaboration VIII 2013). This allows for a 
one-to-one comparison of the simulated and data maps. Each 
simulation leads to the generation of per-detector and per- 
detector-set maps for each survey, and for the nominal mission 
on which the results in this release are based. 



Appendix B: Comparison with a very severe 
selection of data 

This appendix describes the very severe pointing period selec- 
tion (called hereafter superclean) and the impact at the map and 
angular power spectrum levels. 



B.1. Data selection 

Nominal TOI processing is used; only the selection of rings (i.e., 
data obtained during a stable pointing period) is changed and de- 
scribed here. The selected rings include only those showing suf- 
ficiently stationary noise, within strict bounds. To keep enough 
data, the full mission is used for this test, and not only data from 
the nominal 15.5 months. We list here the applied criteria: 

- Normally discarded rings. Of course, rings that are already 
discarded remain so. 

- 4 K cooler lines. Resonant rings are the ones for which a 
4 K cooler line (Sect. 3.6) coincides with a harmonic of the 
spin frequency. The 4 K line correction comes from an inter- 
polation that could introduce some problems. The 10 Hz line 
is the most annoying, so we discard all the 10Hz-4K line 
resonant rings. This criterion is common to all detectors. 

- Baseline jump correction. Rings where a jump correction 
(Sect. 3.8) was performed are discarded. 

- Bursts in the noise. Any anomaly in the noise level (as 
measured by the total noise NET) can be discarded. From 
Sect. 3.10, these are due to either a baseline drift or a small 
jump. The basic threshold used here is ±3 % of the median 
noise per ring. 

- Two-level noise. As described by Planck Collaboration ES 
(2013), the following bolometers exhibit a larger than usual 
noise level in a specific range of rings indicated between 
brackets: 

- 13_353_2 : [4994, end] 

- 23_353_3a : [5811, 9057] and [18890, end] 

- 50_143_3a : [9410, end] 

- 54_353_5b : [2572, end] 

The corresponding ranges are entirely discarded. 

- Integration time per ring. It is possible that the deglitching 
procedure might perform diff'erently for rings with a signi- 
ficantly smaller integration time than average. Hence we dis- 
card rings with an integration time less than 35 minutes. This 
criterion is common to all detectors. 

- Planets. Rings including a big planet (i.e., flagged as in 
Sect. 3.1) can perturb the mapmaking solution if long tails 
are unseen. This criterion is common for detectors in the 
same row of the focal plane. 

- Medium and small RTS. All rings known to be afl'ected by 
RTS (Sect. 3.10) are discarded. 

- Fifth survey. Survey five is fully discarded because of solar 
flares, unusual integration time and end-of-life operations. 
This criterion is common to all detectors. 

- PSB. A logical-OR of both discarded ring lists of the two 
bolometers of a PSB pair is performed. 

B.2. Resulting ring statistics 

Each criterion has a difl'erent impact on both global statistics and 
those of individual bolometers. 

A logical-OR of all of the criteria is shown in Fig. B.l. Black 
pixels indicate that the corresponding rings of the correspond- 
ing bolometer are not used for the superclean subset, whereas 
the red ones are those that were already discarded with the old 
criteria. Several features are notable, such as difl'erences between 
the highest and lowest frequency channels, and the fact that some 
criteria afl'ect all the bolometers. 

The global impact of this ring selection is that about 41 % of 
the rings are discarded for the full mission: about 35 % for the 
CMB channels, and 45 % for the Galactic ones. With the normal 
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Figure A.l. Top row: data difference maps, for the 143- la bolometer. Bottom row: differences from a simulation including the effect 
of the ADC non-linearity and 4K lines. Left column: S1-S2 survey difference for maps made without destriping. Right column: 
S1-S2 differences after the destriping used to create the delivered data. The match is sufficiently good that one can use such 
simulations to estimate the combined effect of the non-linearity of the ADC on the delivered data. 



criteria, only about 1 % of the rings are discarded for the full 
mission (0.2 % in the nominal mission). 

In order to distinguish the impact of a simple decrease of 
the observation time and the impact of the use of the superclean 
selection, we check the difference between discarding 35 % of 
the rings with the criteria above, and discarding 35 % of the rings 
randomly for the 143 GHz channel. 



B.3. Impact of superclean ring selection on Map making 



We use the standard pipeline to create maps for each frequency. 
Along with the global reduction of hit number, the superclean 
ring selection will create stripes on the hit-maps, due to the tem- 
poral structure of some criteria. The sanity check with the ran- 
domly selected rings will have only the impact of the reduc- 
tion of the hit number. This can be seen on the right column 
of Fig. B.2: the number of hits in the superclean maps (middle) 
are not as homogeneous as in the normal map (top) or the ran- 
dom map (bottom). On the left panel, the top map represents the 
masked intensity of the normal map, whereas the middle and the 
bottom maps correspond to the difference with the superclean 
and the random map, respectively. (The differences have been 
whitened with 1/nutcs = 1 /^hitNormai + 1/^hitsevere. where ^hit is 
the number of hits of a map.) 



B.4. Impact of superclean ring selection on spectra 

We then compare the spectra for the superclean ring se- 
lection and the normal ring selection. We use cross-spectra 
of offset-corrected half-ring maps, FIRSTxLAST. The off- 
set correction is done in order to avoid discrepancies at 
low-^ due to mask leakage. To avoid contamination from 
the galaxy, the masks applied to compute these offsets 
depend on the frequency: f,^y= [70,70,40,40,20,20]% for 
[100,143,217,353,545,857] GHz. Notice that, as the same sky 
mask is used with the normal and superclean maps, the errors 
bars do not include cosmic variance. 

Actually, in the Gaussian case, the expected variance of the 
difference is induced by noise, and can be computed as 



■^expected " Q j " f^^^ (2{ + I) ' 

where is the pseudo- spectrum for the superclean selection, 
is the pseudo- spectrum for the normal selection, A^^ is the 
noise pseudo- spectrum of the normal map, and g = 0.65 (or 0.55 
for Galactic channels) is the average ratio of the hit numbers in 
the superclean map with the hit numbers in the normal map. This 
formula holds if the noise of the FIRST and LAST maps have 
the same pseudo- spectrum and are totally uncorrected, and if 
the only difference between the normal and the severe selection 
is the number of hits in them. 
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Figure B.l. Example of normally discarded rings (red) and severely discarded rings (black) for 201 consecutive pointing periods 
(the abscissa has an arbitrary offset). 



Figure B.3 shows those normalized differences for all of the 
channels. The difference is plotted in units of crexpected(C^ - C^)^ 
so that the expected error bars are ± 1 . This is the case from 
^ = 500 to ^ = 2500 for all frequencies. At 143 GHz and 
217 GHz, there is a rise in the pseudo- spectrum for i < 500. 
Such a rise is not seen for the pseudo-spectrum computed from 
the random discarding of 35 % of the rings. However, such a 
feature is far below the cosmic variance (at the percent level). 
We conclude that the level of non-stationarity of the noise in the 
nominal mission data is negligible and does not impact the cos- 
mological results. 
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Figure B.2. Top-left: 143 GHz intensity. Middle-left: difference with superclean. Bottom-left: difference with random selection. 
Those differences have been whitened with the corresponding hit number, shown in the right column. For the superclean selection, 
striping of the hit numbers is more visible.. 
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Figure B.3. For each of the HFI channels, the difference of the normal spectrum and the spectrum obtained with a more severe 
selection of rings, is shown in units of the expected standard deviation. If our assumptions are correct, the residual noise is supposed 
to have a Gaussian distribution centred on zero with a dispersion of one. The red points for the 143 GHz correspond to a sanity 
check with a random discarding of rings (the points here are band- averaged in bins of width ~ 31). 
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